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Summary 

Computer subroutines for approximating un- 
steady aerodynamic forces in the time domain and 
for obtaining a set of constant-coefficient first-order 
differential equations are described. Information is 
provided for using and maintaining the subroutines. 
The unsteady aerodynamic forces are approximated 
by means of either a weighted-matrix Fade method 
or a least-squares method. The subroutines are capa- 
ble of applying a gust spectrum to the data prior to 
calculating the approximations. This technique can 
be used to enhance the accuracy of approximations 
to unsteady aerodynamic forces resulting from ran- 
dom excitation. The approximations may be applied 
over a truncated set of modes at a selected number 
of reduced frequencies. Once the approximations to 
the unsteady aerodynamic forces are calculated, the 
coefficients are combined with the structural matri- 
ces defining the dynamic response of the vehicle to 
formulate a set of four matrices. These four matri- 
ces, or the system quadruple, are the coefficients of 
a set of first-order differential equations and output 
equations that represent the dynamic response of the 
vehicle in the time domain. The user is provided with 
a method of defining output equations that are de- 
rived from a linear combination of the system states. 
Examples are presented to demonstrate the use of 
the subroutines. 

1 Introduction 

The computer subroutine package described 
herein, named UNAERO, was written to imple- 
ment and verify the method reported in reference 1. 
The package contains a collection of FORTRAN IV 
subroutines specifically written for use on a GDC 
CYBER computer system. UNAERO is intended to 
aid the engineer in obtaining a set of linear constant- 
coefficient differential equations of motion that repre- 
sent the dynamic behavior of a flexible airplane. This 
document serves to instruct and aid the user of UN- 
AERO in installing the package, generating the in- 
put, understanding the output, and maintaining the 
package. 

Unlike other computer programs, such as that de- 
scribed in reference 2, UNAERO was designed to 
be incorporated into the user’s existing software. 
UNAERO returns to the calling program the coef- 
ficients of the differential equations of motion as a 
set of four matrices, or the system quadruple. The 
user’s program can use the system quadruple to pro- 
duce the desired data, such as the eigenvalues of the 
system matrix or the singular values of the return 
difference matrix. 


The task of generating the equations of mo- 
tion of the vehicle is divided into the following two 
parts; task I, the generation of the coefficients of the 
unsteady-aerodynamics approximation, and task II, 
the generation of the system quadruple. Within 
UNAERO, either task may be performed indepen- 
dently of the other. The standard way of using UN- 
AERO is to perform task I and then to use those re- 
sults in task II. However, externally generated task I 
data may be used in UNAERO for task II. Within 
task I, the user has the option of using either a 
least-squares approximation method or the weighted- 
matrix Pad4 method. In either method the user can 
weight the error between the approximations and the 
unsteady-aerodynamic-force data to achieve the de- 
sired accuracy of fit at a specified frequency. In ad- 
dition, the aerodynamic-force data representing ran- 
dom inputs can be modified by applying a gust spec- 
trum before the approximations are calculated. UN- 
AERO has a feature which allows the user to reduce 
the number of structural modes or the data or both 
before the approximations are calculated. Within 
task II, sensor equations can be added to the system 
equations through the use of subroutines supplied by 
the user. 

The following is a brief outline of the major sec- 
tions of this paper. A description of the methodol- 
ogy of the program functions along with the options 
and guidelines is presented in section 3. Section 4 
is a technical description of the equations that form 
the basis of the subroutines within UNAERO. In sec- 
tion 5, the command input file for the execution of 
the program is described. To define the vehicle be- 
ing analyzed, UNAERO uses two data bases, one for 
task I and one for task II. These data bases are de- 
scribed in section 6. Program printed output and 
associated options are described in section 7. Sec- 
tion 8, a large section concerned with the structure 
of the program, follows section 7. The majority of 
this material need only be referred to by those read- 
ers who are installing and maintaining the program. 
Finally, examples of the program use are presented 
in section 9. 

If the reader is using this document for instruc- 
tion on the use of the program, sections 3, 4, 5, 
7, and 9 should be read in their entirety, section 6 
should be studied up to and including the description 
of the discrete-frequency data base, and the part of 
section 8 on program use should be read. The in- 
formation required to bypass the aerodynamic-force 
approximation calculations and to formulate a set of 
system equations from data supplied by the user’s 
own technique is contained in the description of the 
continuous-frequency data base in section 6, in the 
appropriate information from sections 4, 5, and 7, 


and in the last part of example 1 in section 9. The 
information presented in appendix A describes the 
interface subroutines PRESEN and SENSOR for use 
in adding output equations to the generated system 
quadruple. Subroutine USER is described in ap- 
pendix B. This subroutine is used to modify the value 
of the aerodynamic force before the approximation is 
calculated. 

2 Symbols 


A 

matrix of coefficients, a function of 
the reduced frequency ki used in 
solution of aerodynamic-force 
approximation (see eq. (12)) 

At 

pseudoinverse of A 

A 

weighted A 

A0,A1,A2 

coefficient matrices of aerodynamic- 
force approximation of equation (6) 

b 

vector of constant coefficients 
used in solution of aerodynamic 
approximation (see eq. (13)) 

b 

weighted b (see eq. (15)) 


mth element of b 

C 

reference length 

D 

generalized damping matrix 

Dm 

coefficient matrices of aerodynamic- 
force approximation of equation (6) 

I 

identity matrix 

•^j 

objective function (see eq. (9)) 

j 


K 

= (n)(iV + 2) 

K 

generalized stiffness matrix 

k 

reduced frequency, 0Jc/2V 

k 

specific reduced frequency 

L 

gust characteristic length 

M 

= {n){Nk) 

M 

generalized mass matrix 

N 

order of the approximation 

Nk 

number of reduced-frequency values 

n 

number of modes used to calculate 
the approximation 

P 

an element of P 


P,P(yo;) 

matrix of numerator polynomi- 
als of unsteady-aerodynamic- 
force approximation 

Pijiik) 

fth,yth element of P at the reduced 
frequency ki 

P m 

matrix of coefficients of P(joj) (see 
eq. (5)) 

Q{i^) 

unsteady-aerodynamic-force matrix 

Q[j^) 

approximation to the unsteady- 
aerodynamic-force matrix 


ith,yth element of Q{juj) at the 
reduced frequency ki 

QijUk) 

zth,yth element of Q{juj) evaluated 
at reduced frequency ki 

q(iw) 

vector of generalized coordinate 

R, R(yo;) 

matrix of denominator polynomials 
of unsteady-aerodynamic-force 
approximation 

R(iw) 

matrix R(yu;) with leading coeffi- 
cient equal to I 

RjUk) 

yth element of R(yw) evaluated at 
reduced frequency fcj 

Rj 

yth element of matrix R 

R, 

matrix of coefficients of R(yu;) (see 
eq. (4)) 

r 

coefficient of the polynomial R 

s 

semispan 

u 

control state vector (see eq. (27)) 

V 

reference velocity 

w 

weight matrix defined by equa- 
tion (14) 

w 

real diagonal matrix (see eq. (14)) 

Wm 

an element of W 

X 

system state vector (see eq. (25)) 

X 

solution vector used in calculation 
of unsteady-aerodynamic-force 
approximation (see eq. (13)) 

Xij 

an element of x 

y 

system output vector (see eq. (25)) 

z 

error between the approximation 
and the data value (see eq. (17)) 
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z complex conjugate of z 

2 pV“^ dynamic pressure of the free stream 

a root-mean-square gust intensity 

bj frequency of oscillation 

II II norm, square root of the scalar 

product 

Subscripts: 

af state vector used to model 

unsteady-aerodynamic-force 
approximation 

C control-surface mode 

G gust mode 

m state vector used for generalized 

coordinate 

5 structural mode 

W weighted norm 

6 S, G, or C 

3 Discussion of the Programmed Options 
and Guidelines 

This section is concerned with the description 
and rationale of the different options that can be 
invoked within UNAERO. As mentioned previously, 
UNAERO consists of two major parts, or tasks. The 
first task (task I) is to calculate the parameters 
of the unsteady-aerodynamic-force approximations, 
and the second task (task II) is to formulate the 
equations of motion of the vehicle from these cal- 
culated parameters. In performing these two tasks, 
there are several options that are available. Some of 
the options aid the user in a search for the smallest 
set of dynamic equations that adequately represents 
the aircraft’s dynamic response. Other options are 
for the efficient use of computer resources. The op- 
tions that control the calculation of the parameters in 
the approximations (task I) are described first. The 
options that control the formulation of the dynamic 
equations of motion (task II) are described subse- 
quently. 

3.1 Task I: Options and Guidelines 

In the development of the methodology for cal- 
culating the parameters of the approximations to 
the unsteady aerodynamic forces, the assumption is 
usually made that the dynamics of the aircraft can 
be represented accurately by a finite number of air- 
craft modes. Furthermore, the generalized forces 


acting on the aircraft are a linear combination of 
forces arising from the motion of rigid-body or flex- 
ible modes, motion of control-surface modes, and 
motion of forces resulting from an external random 
disturbance. These three sources are hereinafter re- 
ferred to as the structural modes, control modes, and 
gust modes, respectively. For the structural modes, 
UNAERO requires the generalized mass, damping, 
and frequency in the data base. 

UNAERO makes no distinction between the forms 
of the unsteady-aerodynamic-force approximations 
used for any of the structural or control modes. The 
form of the approximation that is used for the struc- 
tural and control modes is the following rational 
transfer function: 

E Ukfai 

1=0 

E 

1=0 

with the order of the numerator polynomial K be- 
ing two greater than the order of the denomina- 
tor polynomial N. The additional derivatives that 
are required for the structural mode approximations 
are obtained from the equations used to represent 
the motion of the generalized modes. These equa- 
tions are generated with the generalized mass, damp- 
ing, and frequency data. For the control modes, 
UNAERO expands the number of control inputs to 
be three times the original number of control modes. 
The additional inputs represent the first and second 
time derivatives of each control deflection. The form 
of the aerodynamic-force approximation used for the 
gust mode is identical to the form used for the struc- 
tural mode except that it is restricted to be a causal 
function. This form of the approximation for a gust 
mode does not require additional inputs that repre- 
sent the derivatives of the gust mode (i.e., the order 
of the numerator K does not exceed the order of the 
denominator polynomial N for the approximation). 

For all unsteady-aerodynamic-force approxima- 
tions that are dependent on the dynamics of a mode 
(i.e., all elements of a column of the generalized force 
matrix), the user has the option of using the least- 
squares method or the matrix Fade method to solve 
for the denominator coefficients of the approxima- 
tion. The recommended procedure for calculating 
the approximation is to use the Pad4 method first 
and to adjust the parameters of the approximation 
(order of the approximation, number of modes, etc.) 
until the desired degree of accuracy is obtained. Then 
the results of the Fade method are used as the values 
of the denominator coefficients for the least-squares 
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method. This procedure can result in a substantial 
savings in the resources used for an analysis of a 
vehicle. The savings are obtained because once de- 
nominator coefficients are determined with the Fade 
method, they will be constant. Since these coeffi- 
cients are constant and known, they can be used 
as input to the least-squares method to produce ap- 
proximation functions with identical numerator coef- 
ficients but which require less resources to calculate. 

A special case has been added to the least-squares 
method to allow a zero-order approximation to the 
aerodynamic forces. An approximation of zero order 
has no lag terms, and therefore no additional states 
are required to approximate the unsteady aerody- 
namic forces. In this case, the number of system 
states in the resulting equations of motion is two 
times the number of modes in the structural model. 
Although the approximations of zero order generate 
functions that produce large errors, the error at any 
value of reduced frequency can be minimized through 
the use of a weighted least-squares solution for the 
coefficients of the approximations. An example in 
section 9 illustrates this option. This technique of 
using weights to selectively reduce the error of the 
function can also be used with the Pad^ method and 
the least-squares method when the order of the ap- 
proximation is greater than zero. 

The data being used to generate the approxima- 
tion functions can be adjusted so as to be a subset 
of the original data. If this option is used, a num- 
ber of systems can be generated and evaluated as to 
which modes make a significant contribution to the 
dynamics of the problem being investigated. Also, 
the number of reduced frequencies where data are 
used to calculate the approximating functions can be 
reduced so that the resulting function is based on a 
frequency range that is significant for the mode be- 
ing approximated. This could also be accomplished 
with a value of zero for the weight at the desired 
value of reduced frequency. But it would not be as 
efficient as removing the data value from the data 
base used to calculate the approximating function. 
Reducing the number of modes used to represent the 
dynamics of the vehicle or the number of data val- 
ues used or both will reduce the amount of computer 
resources required to calculate the coefficients of the 
approximation. 

Based on experience gained in using UN AERO, 
it has been observed that approximating gust modes 
accurately is a more difficult task than approximat- 
ing other types of modes. In an effort to improve 
the quality of the gust mode representation, there 
is an option in UNAERO which allows the modifi- 
cation of the aerodynamic data (by a gust attenua- 
tion function) before the approximation coefficients 


are calculated. The attenuation can be based on the 
magnitude of a Von Karman or a Dryden gust spec- 
trum or represented by a FORTRAN subroutine sup- 
plied by the user. The details needed for writing this 
subroutine are contained in appendix B. Because of 
the method by which this option is implemented, the 
data being modified are not limited to the gust mode 
data and can be used for any mode. 

3.2 Task II: Options and Guidelines 

When UNAERO is being used to formulate the 
equations of motion from data that have been previ- 
ously approximated, there are two options that are 
available to the user. In the first option the refer- 
ence condition at which the equations of motion are 
formulated must be specified. This reference condi- 
tion, or set point, is defined and must be redefined 
each time the equations of motion are formulated. 
Of the many different ways that the reference con- 
dition could be specified, for UNAERO a charac- 
teristic length, velocity, and dynamic pressure are 
required. This set of parameters allows the equa- 
tions of motion to be conveniently reformulated at 
another dynamic pressure. One of the first checks 
on the formulation of the equations of motion can be 
an inspection of the system matrix with the dynamic 
pressure set to zero. This ensures that the general- 
ized mass, damping, and stiffness are entered in the 
data base correctly. 

In the second option the user may change the 
coefficients of the output equations, thereby allowing, 
for example, an analysis of the effect of changing the 
location of a sensor. Since adjusting the coefficients 
of the output equations requires the interpolation 
of the modal displacements, the user must supply 
the subroutines PRESEN and SENSOR to use this 
option. The details of these two subroutines are given 
in appendix A. 

4 Description of Programmed Equations 

The equations within UNAERO are described in 
this section. This description is a supplement to the 
technical description of the methods in reference 1 
and is not intended to provide the detailed techni- 
cal development that was presented in that refer- 
ence. The equations of motion that describe the air- 
plane, written in the frequency domain, are described 
first. The imsteady-aerodynamic-force approxima- 
tion functions are then described, followed by a de- 
scription of the objective function and the equations 
used in the optimization procedures to minimize that 
objective function. The time domain equations that 
represent the dynamic response of the aircraft are 
then presented. 
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4.1 Task I Equations 

In the development of the equations of motion 
for an elastic airplane, the motion is assumed to be 
a perturbation from straight and level flight. The 
equations of motion in the frequency domain, which 
were developed in equations (1) to (10) of reference 1, 
are summarized as follows: 


{jujfM + {juj)D + K + -pV^Qiju) 


q(jo;) - 0 


( 1 ) 

where M, D, and K are the generalized mass, gener- 
alized damping, and generalized stiffness matrices, 
is the dynamic pressure of the free stream, 
Q{juj) is the unsteady-aerodynamic-force matrix, 
and q(yo;) is the vector of generalized coordinates. 
To develop the equations that are used in UNAERO, 
equation (1) is rewritten with the aerodynamic forces 
separated to take into account the notation for the 
different modes that were described earlier. The 
modified equation (1) is 


which best fits the jth column of the matrix Q{ju)) 
at a number of values of reduced frequency. Each el- 
ement of the matrix P(j w) is a polynomial of degree 
N + 2 in joJ. Each diagonal element of the matrix 
R(ya;) is a polynomial of degree N in ju. Because 
R(yti;) is a diagonal matrix, each column of the ap- 
proximation matrix Q(jw) will have the same denom- 
inator polynomial. Equation (4) represents the de- 
nominator polynomials and equation (5) represents 
the numerator polynomials of equation (3) rewritten 
so that all the numerator coefficients appear as con- 
stant matrices. 

Equations (4) and (5) are the internal representa- 
tion of the approximation for the UNAERO package. 
The elements of the denominator polynomial R(yo;) 
are defined in table 1. The polynomial used is differ- 
ent for each order of the approximation because of 
the technique used to solve for the elements of the 
matrix R(jo;). The form of equation (5) that ap- 
pears the most in the literature and that is used in 
UNAERO during the task II phase when formulating 
the equations of motion is written as 


[(jw)^ M + (jw)D + K + ] qsO'w) 

“ “ (2) 
The subscripts S, C, and G denote a general- 
ized force or coordinate associated with a structural 
mode, a control mode, and a gust mode, respectively. 

The equation used to represent the approxima- 
tions of order N to the unsteady-aerodynamic forces 
due to motion of the structural modes or the control 
modes can be written as equation (3). 


Q^O'w) = AO^ + (jw) Al{ + {jwf A2g 


N-1 


~ / r \ m-N 


,m=0 


R«0w)-i(6) 


where the matrix R(yw) is defined as 

^ \ -N 




(7) 


where 6 represents either S or C, 


(3) For forces resulting from the random or gust excita- 
tion, the approximation function is written as 


N A-1 _ I 

Rs(ji^) = 1 + (^) ( 4 ) 

1=0 


m=0 

The matrix Q{ju>) is the approximation in the fre- 
quency domain of the unsteady aerodynamic forces, 
P{juj) is a matrix such that the fth,yth element rep- 
resents the numerator polynomials of the unsteady- 
aerodynamic-force approximation, and R(ju;) is a di- 
agonal matrix such that the yth diagonal element 
is the denominator polynomial of the approximation 



Qg(J^) = A0(j 

fA-l 

+ 


/ r 


Lm=0 




( 8 ) 


The matrices AO, Al, A2, and Dm in equations 6 to 
8 and the matrix R/ in equation (4) are the matrices 
computed by UNAERO. Note that equation (8) is 
equivalent in form to equation (6) with matrices Al 
and A2 set to zero. 


4.1.1 Least-squares approximation equations. The 

objective function, or the function that is minimized 
by calculating the coefficients of the approximations 
with either the least-squares method or the Fade 
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method for the jth column of the unsteady aerody- 
namic forces, is: 

1 n Nk 

(9) 

^ i=i 1=1 

where n is the number of modes used in the ap- 
proximation, Nk is the number of reduced-frequency 
values used to calculate the approximation, Qij{jki) 
is the value of the ith,yth generalized aerodynamic 
force at the reduced frequency fc/, and Qijijki) is 
the value of the ith,jth approximation function eval- 
uated at reduced frequency fc;. 

For the least-squares approximation method, 
equation (9) is minimized by solving for the numer- 
ator coefficients with the method described in ref- 
erence 1 and modified so that weights on the error 
between the approximation and the calculated force 
at specified frequencies can be applied. Note that 
when the denominator coefficients of the approxima- 
tion are constant, which is the case for the least- 
squares method, then the solution to the minimiza- 
tion of the following equation will also be the solution 
that minimizes equation (9): 


can be rewritten as 

J, = i||b-Ax||2 (13) 

Note that the summation of the scalar norm of 
equation (10) is equivalent to the vector norm of 
the algebraic quantity in equation (13). The real 
vector X contains all the coefficients of the numerator 
polynomials that are the solution of the least- 
squares optimization procedure. 

To modify equation (13) for the weighted opti- 
mization procedure, define a weight matrix as 

W = W'W (14) 

where W is a real M xM diagonal matrix of weights 
and the prime indicates the transpose of the matrix. 
Then, by defining 


b = Wb 1 
A = WA J 


equation (13) becomes 


(15) 


i n Nk 

= 2 E E \\QijUkl)Rjijki) - Pij{jki)f (10) 

i=l 1=1 

where P^j is the numerator polynomial for the ith,jth 
element and Rj is the denominator polynomial for 
the jth column of the approximation; both are eval- 
uated at the reduced frequency of ki- If the elements 
of the summation over the indices / and i of equa- 
tion (10) are grouped together to form a vector b of 
size M = {n){Nk), the mth element of b would be 
defined as 

bm = Qij{jki)Rj{jki) m = n{l-l) + i (11) 

The coefficients of the numerator polynomials for 
the jth column are also grouped into a vector x, 
of length K = {n){N + 2), such that when this 
vector is premultiplied by an M x K matrix A, 
the mth element of the resulting vector is the fth 
numerator polynomial evaluated at the ki value of 
reduced frequency, or 

Pijijkl) — (Ax)^^^jj element 

The matrix A is a complex matrix formed from the 
discrete values of reduced frequency being used to 
calculate the approximation. Thus, equation (10) 


Jy = i||b-Ax||2 = i||b-Axf„ (16) 

The vector x is found such that it minimizes equa- 
tion (16) by the least-squares method described in 
the appendix of reference 1. 

4.1.2 Fade approximation equations. In order to 
minimize the Fade objective function, a nonlinear op- 
timization algorithm is used that requires analytical 
gradients. To obtain these gradients, equation (9) is 
rewritten as 

Jj = iz'Wz (17) 

where z is a complex vector of length M = (n) {N^) 
with an mth element of 

— QijiJkl) ~ Qijijki) (18) 

where the subscript m is i + n{l — l), z is the complex 
conjugate of z, and n is defined in equation (9). 
The mth element of z represents the error between 
the data and the approximation at the Zth value of 
reduced frequency for the ith mode. The term W is a 
diagonal matrix of weights such that the mth element 
is the weight on the error for the fth mode at the fth 
value of reduced frequency, and N^. is the number of 
reduced-frequency values being used to calculate the 
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approximations. Equation (17) can be rewritten as 


1 ^ 

Jj = - ^ + lm{zrnf 


W„ 


(19) 


m=l 


where Re(2m) and Im(z,„) are the real and imaginary 
parts of the mth element of z, Wm is the mth 
element of the matrix W defined in equation (17), 
and M = (n)(7V^). Taking the partial derivative of 
Jj with respect to a parameter r of the denominator 
polynomial yields 


„ M 

17 = E (^) +I™(^»n)Im [^)] Wm 

m=l 

( 20 ) 

IfRy is an element of the matrix R in equation (4) 
(for which the coefficients are defined in table 1) 
and P is the jth column vector of the matrix P 
in equation (5), then the column vector that is the 
partial derivative of z with respect to the parameter r 
can be calculated from a relationship of the form 


dz 

dr 


P dP 





( 21 ) 


where the appropriate elements of the polynomials 
P and Rj must be used to be consistent with the 
element of z for which the derivative is being cal- 
culated. The relationship for dP/dr is obtained 
from the solution of the least-squares minimization 
of equation (13) as 



Equations (20) to (22) and the denominator poly- 
nomials defined in table 1 are used to calculate the 
gradient of equation (17), which is then used in the 
optimization procedure to calculate the denominator 
coefficients of the approximation function. 


4.1.3 Gust attenuation equations. In order to 
modify the calculated unsteady-aerodynamics gust 
data with a gust attenuation function before the 
approximation function is calculated, the following 
equations are evaluated to obtain the attenuation to 
be applied. 

Von Karman attenuation function: 


^vk(^) = 


l-|-(32/3)[1.339(L/c)fc]^ 
{1 + 4[1.339(L/c)A:]2}“/® 


(23) 


Dryden attenuation function: 




1 + 24[{L/c)kf 


(24) 


The quantity a is the root-mean-square gust intensity 
and L is the characteristic length, or scale, of the 
gust excitation. Equations (23) and (24) are based 
on equations (6) and (8), respectively, of reference 3 
and have been modified to account for the differences 
in the Fourier transform used in references 1 and 3. 


4.2 Task II Equations 

The time domain equations of motion, or sys- 
tem quadruple, are formulated in the following man- 
ner: (1) the approximations for the unsteady aero- 
dynamic force (eqs. (3) to (6)) are substituted into 
the frequency domain formulation of the equations 
of motion (eq. (2)), and (2) the resulting equations 
are transformed into the time domain by perform- 
ing an inverse Fourier transformation. (The details 
of those transformations are in ref. 1.) The time do- 
main equations are then in the following form: 


X = AX -I- Bu 1 
y =2 CX -I- Du j 


(25) 


where A, B, C, and D form the system quadruple 
that defines the equations of motion, and X is a 
state vector composed of all the states to represent 
the dynamics of the vehicle and has the following 
definition: 

{ Qm 1 

qm I (26) 

where is a vector containing the generalized 
coordinate of any mode. The vector is the 
time derivative of the vector Qm. The vector Qgf 
is a vector formed by grouping the states used to 
approximate the unsteady aerodynamic forces. 

The elements of the vector Qgf are arranged such 
that the states that are used for an individual ap- 
proximation are in sequence. For example, if the 
approximation of the unsteady aerodynamics due to 
motion of the first mode is third order and the ap- 
proximation for the aerodynamics due to motion of 
the second mode is second order, then the first three 
elements of Qgf will be the state variables used in the 
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approximation of the aerodynamics due to motion of 
the first mode and the next two will be used in the 
approximation for the aerodynamics due to motion of 
the second mode. Additional unsteady-aerodynamic- 
force states due to the motion of modes will be added 
in a similar manner. Note that the order of the states 
is different from the order of the states that is used 
in equation (16) in reference 1. The control vector u 
is expanded to provide the required input as follows: 


' U\ 
ill 
ill 


llm 
ilm 
ilm 
. Ug . 


(27) 


where, for this illustration, there are m aerodynamic 
controls ui to Um and one gust input Ug. The 
output vector y is defined by the user and can be any 
linear combination of the state vector, the derivative 
of the state vector, and the control vector. (See 
appendix A.) 


5 UNAERO Command Input 

There are two distinct sets of program input. The 
first set, the command input, selects the different op- 
tions of the program and is described here. (The 
second set, the data bases, defines the unsteady- 
aerodynamic-force model and is discussed in sec- 
tion 6.) Examples of the command input are given 
in figure 1 and are executed by the program in fig- 
ure 2. The syntax of the command input file was 
designed to be self-documenting in that the char- 
acteristic parameters of the model being generated 
can be understood from examining the input data 
without the need of a reference manual. Also, an 
effort was made to allow the user an easy method 
for changing the mutually exclusive options. For ex- 
ample, to change from a Fade approximation to a 
least-squares approximation, the analysts need only 
change the phrase “FADE” to “LEAST SQUARES” 
on the appropriate data record. 

5.1 Key Elements 

Throughout this discussion, the term “record” 
refers to one card image. Each input record has space 
allocated for a description of the data that are con- 
tained on the record. A position-dependent portion 
of this description will be referred to as the “key el- 
ement.” The key element is used to differentiate the 
command input data by UNAERO. Exceptions to 


this are data records that can have a variable num- 
ber of items, such as data records that define the 
modal selection array or the frequency selection ar- 
ray. For these exceptions, the description and the 
key element are contained on the previous record. 
The key elements must appear in the correct posi- 
tion in the record relative to the start of the record. 
Figure 3 is a listing of all the legal command input 
records. The relative position of each key element 
can be found by comparing the column position of 
the key element with the first and second entries in 
the figure. For example, in figure 3 the key element 
of the first command input, “E LE,” begins at the 
ninth relative position from the start of the record. 
If the key element does not agree with the expected 
key element or if there is an error while decoding the 
data, such as an illegal character in the data field, 
an error message is generated. The key elements of 
figure 1(a) are underlined in figure 3. 

5.2 Syntax Diagrams 

The command input is best described by the use 
of the syntax diagrams of figure 4. The diagrams are 
used to construct the input records for UNAERO. 
Starting at the upper left arrow of each diagram, the 
path through the diagram is indicated by the arrows. 
A path through the diagram defines a syntactically 
correct input into UNAERO. Each box within the 
diagram references by name another diagram, which 
is then used to define the expected input. For 
convenience, the figure number that this box refers 
to is in the lower right-hand corner. For example, the 
word “Block” in figure 4(a) refers to figure 4(b). The 
word “Comment” in figure 4(b) refers to figure 4(c). 
This process continues until a terminal symbol is 
reached. A terminal symbol represents actual data 
input into the program and is indicated in figure 4 by 
a rounded enclosure. For an example, see figure 4(d). 
In figure 4, the key element in each of the terminal 
symbols is underlined. Data formats are enclosed in 
parentheses and must be strictly adhered to. 

5.3 Example of Using Syntax Diagrams 

Consider for this example the following single 
command statement: 

REFERENC E LE NGTH (FIO.O) 

The character sequence “E LE” is the key element 
beginning at the ninth relative position from the 
beginning of the record. Also, there is a data field 
starting one space after “LENGTH” which is read by 
the program with a format of FIO.O. The path for 
this example is as follows: 

1. The command input starts with figure 4(a). 
The word “Block” refers the path to figure 4(b). 
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2. In figure 4(b), the path for this example con- 
tinues to figure 4(d), “Statements for calculating 
approximation.” 

3. In figure 4(d), the path is referred to fig- 
ure 4(f), “Define set point.” 

4. In figure 4(f), the terminal symbol is reached. 

From figure 4(f), the program expects the next state- 
ment to be 

REFERENC E VE LOCITY (FIO.O) 

This input statement defines the reference velocity. 
The units must be consistent with the reference 
system that was used in developing the data base. 
Note that, according to figures 4(b) and 4(c), if the 
original record had been 

USE APPROXIMATION ON VDATA FILE 

instead of “REFERENCE LENGTH,” the “REFER- 
ENCE VELOCITY” statement would have gener- 
ated an error. Examples of syntactically correct in- 
put can be found in figure 1. 

6 UNAERO Data Bases 

This section describes the second set of input 
to UNAERO, the data bases. Since the input de- 
scribed here is considered to be a data base rather 
than command input, as described in section 5, the 
file structure is not self-documenting. Allowance 
is made for comments to describe the blocks of 
data. The data bases contain a description of the 
unsteady-aerodynamic-force model of the aircraft un- 
der investigation. 

There are two data bases, one for each of the basic 
tasks. For task I the data base contains information 
describing the unsteady-aerodynamic-force model at 
a specific number of reduced frequencies, as would 
be supplied by an analytical oscillatory aerodynamic 
program. Henceforth this data base is referred to 
as the discrete-frequency data base, or the discrete 
data base. The second data base, used for task II, 
differs from the first in that the aerodynamics of the 
model have been approximated and the coefficients 
of the approximation have replaced the unsteady- 
aerodynamics data of the discrete-frequency data 
base. Only the selected modes which are specified 
in task I are retained in the discrete-frequency data 
base. This form of the data base will be referred to 
as the continuous-frequency data base, or the contin- 
uous data base. The continuous-frequency data base 
can either be generated internally from the discrete- 
frequency data base by execution of UNAERO or 
externally by some other method. If generated ex- 
ternally, the user is responsible for ensuring that 
the data format and structure conform to what UN- 
AERO requires. This provides a convenient method 


for comparing different unsteady-aerodynamics ap- 
proximation techniques. 

6.1 Description of the Discrete-Frequency 
Data Base 

The discrete-frequency data base has all the in- 
formation necessary to calculate the approximation 
to the unsteady-aerodynamics model and to gener- 
ate the continuous-frequency data base. Figure 5 
contains a partial listing of an example of this data 
base. This data base is large even for a small num- 
ber of modes and reduced frequencies. In order to 
reduce the amount of data presented, only the sig- 
nificant portions of the data base are listed. The 
data base contains an unsteady-aerodynamics model 
for the DAST ARW-1 wind tunnel model at a Mach 
number of 0.9. It will be useful, in the following dis- 
cussion, to refer to figure 5. 

The discrete-frequency data base contains manda- 
tory comment records. The comment records are not 
optional and must be present even if they are blank. 
In figure 5, records containing alphanumeric infor- 
mation are the mandatory comment records and the 
numbers in parentheses on these records correspond 
to the numbered items in the following list. These 
numbers in parentheses have been placed to the left 
of the actual record and do not represent actual data 
in the data base. Each numbered item (1) to (8) 
describes a section of the data base. 

(1) The first record of the data base is a comment 
record used to describe the contents of the data 
base. The second record describes the format for 
reading all real-type data arrays in the rest of 
the data base. This record must be a syntac- 
tically correct FORTRAN FORMAT statement 
and must be enclosed in parentheses. In the fol- 
lowing sections, this format is stored in the array 
called IFORM. 

(2) The second record beginning this section de- 
scribes the system parameters that are on the 
next record. The next record following the com- 
ment record contains four integer-type data ele- 
ments in the format of 415. These four system 
parameters are the number of structural modes 
NM, the number of control modes NC, the num- 
ber of gust modes NG, and the number of reduced 
frequencies NRF for which there are aerodynamic 
data. 

(3) The third section of the data base contains the 
generalized mass matrix GM. This is a square ma- 
trix and UNAERO uses the following FORTRAN 
statement for reading GM: 

READ(JSl,IFORM)((GM(I,J),J=l,NM),I=l,NM) 
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(4) The generalized damping coefficients are specified 
in this section. UNAERO reads NM data ele- 
ments with the IFORM format specification. The 
data elements are placed along the diagonal of the 
square generalized damping matrix. 

(5) The natural frequencies (in hertz) for each of the 
structural modes are read in this section. UN- 
AERO reads NM data elements with the IFORM 
format and creates a square generalized frequency 
matrix. 

(6) The reduced frequencies for which the oscillatory 
aerodynamic data were calculated are read in this 
section. UNAERO reads NRF data elements with 
the IFORM format. 

(7) This section contains the real parts of the gener- 
alized aerodynamic forces for the first frequency. 
UNAERO reads these data with the following 
statement: 

READ(JS1,IF0RM)((DUM(I,J),I=1,M1), 

J=1,M2) 

where DUM is the variable array where the data 
are initially stored, Ml is the sum of the number 
of structural modes and the number of controls 
NM-I-NC, and M2 is the sum of the number of 
structural modes, the number of controls, and the 
number of gust modes NM-t-NC-l-NG. 

(8) This section contains the imaginary parts of the 
generalized aerodynamic forces for the first fre- 
quency. The data are read with a statement iden- 
tical to that in item (7). 

Items (7) and (8) are repeated for each reduced 
frequency that is contained in the data base. Op- 
tional data records can be appended to the end of 
the data base for use by the user-defined subroutine 
PRESEN. For a description of the method for adding 
sensors to the model definition, see appendix A. 

6.2 Description of the Continuous-Frequency 
Data Base 

The continuous-frequency data base consists of 
the calculated parameters for the unsteady- 
aerodynamic-force approximation plus all other data 
necessary to construct a system of first-order constant- 
coefficient differential equations at a given flight con- 
dition. Figure 6 contains an example of this data 
base. This data base was generated from the data 
base described in figure 5, but retains only the first, 
second, fourth, fifth, and sixth modes. The gener- 
alized aerodynamic forces were approximated with a 
first-order Fade approximation. As in the description 
of the discrete-frequency data base, only the signifi- 
cant portions of the data base are listed in order to 
conserve space. Each numbered element of the fol- 
lowing list describes a section of the data base that 


has the number enclosed in parentheses to the left 
of the actual data. The comment on the description 
record is not required, but the record is not optional 
and must be present even if it is blank. 

(1) The first record of the data base is a com- 
ment record that describes the contents of the 
data base. (If this data base was generated by 
UNAERO, the first record is the same as the first 
record for the discrete data base.) 

(2) The second section of the data base begins with 
a comment record followed by a data record con- 
taining, in 315 format, the number of structural 
modes NM, the number of control modes NC, and 
the number of gust modes NG that are contained 
in the data base. The next record in this sec- 
tion must be a syntactically correct FORTRAN 
FORMAT statement that describes the format 
for reading the successive data. If the data 
base is generated by UNAERO, this format is 
“(6G12.5).” This format is stored in the variable 
IFORM. 

(3) This comment record indicates that the coeffi- 
cients of an aerodynamic approximation are to 
follow. Since the input data records are the same 
for a structural mode or for a control mode, both 
sets of input are described here. UNAERO uses 
the comment “MODE” to indicate that an ap- 
proximation to a structural mode is to follow and 
the comment “CONTROL” to indicate data for a 
control mode is to follow. The data record follow- 
ing this comment record contains the order of the 
aerodynamics approximation NP in an 15 format. 
UNAERO reads the coefficients of the approxi- 
mation from the next data record or records with 
the following statements: 

IF(NP.NE.O)READ(JD,IFORM)(R(I),I=l,NP) 

READ(JD,IFORM)((P(J,K),K=l,NPP3), 

J=1,NM) 

where NP and NM are defined above and NPP3 
= NP -I- 3. If the order of the approximation 
is zero (NP=0), the R array is not present in 
the data base. The polynomial defined with the 
R coefficients of the R array will be different 
depending on the order of the approximation NP. 
These polynomials are defined in table 1. The 
coefficients of the P array are used to form NM 
polynomials of order NP + 2 with the coefficients 
of the polynomial being assigned in decreasing 
powers of reduced frequency. 

(4) This section defines the approximation for the 
gust. Following a comment record, the first data 
record contains the number of coefficients in the 
denominator NPl and in the numerator NP2 in a 
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format of 215. The quantity NPl is the order of 
the denominator polynomial, whereas NP2 is the 
order of the numerator polynomial plus one. If 
NPl is not zero, the coefficients of the polynomial 
are read by 

READ(JD,IFORM)(R(I),I=l,NPl) 

The coefficients of the denominator, in descending 
order, follow and are read by 
READ(JD,IFORM)((P(J,K),K=l,NP2), 
J=1,NM) 

For the example in figure 6, the unsteady- 
aerodynamics effects of the random gust were ap- 
proximated by a constant. Therefore the values 
of the corresponding R vector are not present and 
there are five elements of the P matrix. (See 
eq. (4).) 

(5) This section defines the generalized mass data. 
Following a comment record, UNAERO reads the 
generalized mass matrix with 
READ(JD,IFORM)((GM(I,J),J=l,NM), 

I=1,NM) 

(6) This section describes the generalized damping 
vector. Following a comment record, UNAERO 
reads NM damping coefficients from the data base 
with the format statement stored in the variable 
IFORM. 

(7) The natural frequencies (in hertz) for each of the 
structural modes are the final items. Following a 
comment record, UNAERO reads NM data ele- 
ments according to the format statement stored 
in the variable IFORM. 

Additional data records can be appended to the 
data base as required by the user-defined subroutine 
SENSOR. If the data base is generated by UNAERO. 
rather than from an external source, these additional 
data records are read by the user-defined subroutine 
PRESEN. See appendix A for an example of using 
the subroutines PRESEN and SENSOR. 

7 UNAERO Output 

UNAERO prints selected output and returns the 
system quadruple to the calling program through the 
calling sequence. The user can direct UNAERO to 
list the command file, the calculated parameters of 
the aerodynamics approximation, the original and 
the estimate of the data, the error of the approxi- 
mation, the value of the objective function, the co- 
efficients of the system equations at the specified set 
point, and diagnostic information. 

Figure 7 contains a sample of the printed output 
from UNAERO with IPRINT set equal to 4 by the 
calling program. Table 2 contains the relationship 


between the value of IPRINT and the resulting out- 
put. In the following list, each numbered item in the 
list corresponds to the circled number in figure 7. (If 
IPRINT had been set to 2 rather than to 4, then 
item 1 would be the only item printed.) 

1. The command input file. 

2. The calculated coefficients of the aerodynamic 
approximation. For the results printed in this 
example the approximation method is the Pade 
method. Therefore the P and R matrices are 
printed. If a least-squares method is used, then 
the R matrix is suppressed. 

3. Indicates that a comparison between the origi- 
nal data and the calculated approximations for 
the indicated mode follows. The first number 
is the mode number in the continuous-frequency 
data base at which the calculated data are stored. 
The number in parentheses is the mode in the 
discrete-frequency data base from which the coef- 
ficients of the approximation are obtained, and it 
corresponds to a column of the generalized force 
matrix. 

4. Indicates that the following information, items 5 
to 10, are the data for the approximation of the 
listed mode number in the continuous-frequency 
data base. This corresponds to a row element of 
the generalized force matrix. 

5. The frequencies at which calculations (listed to 
the right) were made. 

6. The data from the discrete data base that were 
used in the approximation. 

7. The calculated approximations at the indicated 
value of reduced frequency. 

8. The absolute error between the original data and 
the approximations. 

9. The relative error between the original data and 
the calculated approximations, expressed as a 
percentage of the original data. 

10. The square root of the weight that was applied to 
the squared error in the objective function. 

11. Items 4 to 10 represent the printed analysis of the 
approximations to the generalized aerodynamic 
forces that are due to motion of the mode indi- 
cated in item 3. Items 4 to 10 are repeated for 
each mode in the system being generated that the 
generalized force acts upon. In this example there 
are two sets printed, one for each of the modes in 
the system being generated. 

12. The maximum absolute and relative errors for the 
approximation listed at the end of the analysis 
along with the value of the objective function, 
identified as “COST” in the sample printout. The 
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lower the objective function, the more accurate 
the fit. 


itemizes the information contained on the microfiche 
pages. 


13. The system matrices at the requested set point 

printed at the end of the comparison of the data 

and the approximation. 

If this example were complete, UNAERO would 
repeat items 2 to 12 for each approximation function 
being generated. In this example there are three ap- 
proximation denominator polynomials, one for each 
mode of the model and one for the control force. 
Only the data from the first approximation appear 
in figure 7. 

For the convenience of the user, the printed out- 
put has been grouped into three sets chosen to co- 
incide with the three intended uses of UNAERO: 
(1) initial calculation of the aerodynamics approxi- 
mation, (2) model generation with or without printed 
results, and (3) program installation and checkout. 
When UNAERO is used to calculate an unsteady- 
aerodynamics model, the accuracy of the approxima- 
tion needs to be monitored while the parameters of 
the approximation are adjusted accordingly. In this 
application, output set 1 should be chosen. After the 
user is satisfied with the accuracy of the approxima- 
tion, the printed output of UNAERO can be reduced 
to a set that consists only of the command input 
file. This use of UNAERO can arise when systems 
are generated by varying parameters of the equations 
that do not affect the approximation accuracy, such 
as the dynamic pressure or the position of a sensor. 
An option of suppressing all the output exists, such 
as when UNAERO is an intermediate step of another 
algorithm in which the output of UNAERO would 
not be appropriate. The final output option (set 3) 
consists of two print levels that are designed for pro- 
gram installation and checkout. The choice of these 
options requires the user to be familiar with the pro- 
gram modules in order for the printed information to 
be of value. Therefore, set 3 should only be used by 
a user that is installing UNAERO and needs to check 
an intermediate calculation or result or both. The set 
to be printed is controlled by the value of the variable 
IPRINT and is written to the file whose input/output 
unit designator is contained in the variable IW. The 
file that is generated on this unit is intended to be 
routed to a printer with the capability of printing 
records that have a length of at least 130 characters. 

A complete listing with IPRINT equal to 7 is 
present in the microfiche available from the au- 
thor. The information on the microfiche is written in 
columns, with 15 pages per column. The first page 
of the first column is labeled AOl and the last page 
of the first column is labeled OOl. The first page of 
the second column is labeled A02. The following list 


Page 

AOl-BOl 

COl-FOl 

GOl 

HOl-JOl 

K01-N02 

002 

A03 

B03-C03 

D03 

E03-F03 

G03 

H03-M03 

N03 

003 


FORTRAN control program 

CYBER loader map 

Input data that generated this mi- 
crofiche output (also item 1 in fig. 7) 

UNAERO output information for a 
first-order Fade approximation 

CONMIN output 

UNAERO output discussed for item 2 
in figure 7 

UNAERO output discussed for items 3 
to 12 in figure 7 

UNAERO output information for a 
first-order least-squares approximation 
(also item 2 in fig. 7) 

UNAERO output discussed for items 3 
to 12 in figure 7 

UNAERO output information for a 
zero-order least-squares approximation 
(also item 2 of fig. 7) 

UNAERO output discussed for items 3 
to 12 in figure 7 

UNAERO output for task II 

Discrete-frequency data base 

Continuous-frequency data base 


8 UNAERO Structure and Execution 

In this section the interrelation and calling proce- 
dures for the subroutines used within UNAERO are 
discussed. For those users who only require the infor- 
mation necessary to use UNAERO, this information 
is presented first. For those users who will install 
and maintain UNAERO, each major subroutine of 
the software package is discussed in section 8.2. For 
each subroutine a brief description of the function or 
subroutine is given, the calling sequence is defined, 
and the argument list is described. 

8.1 UNAERO Execution 

Since UNAERO is a package of FORTRAN sub- 
routines, the user must provide a program that ini- 
tializes the input/output routines, initializes the vari- 
ables of the labeled common block INOUT, and calls 
the subroutine UNAERO. The program listed in fig- 
ure 2 is an example of a calling program which uses 
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the UNAERO package of subroutines. The subrou- 
tine UN AERO will expect to read the command in- 
put file and data bases from the unit designators 
stored in the variables IR and JD of the labeled com- 
mon block INOUT. On return the system matrices 
that are in the calling sequence will contain the de- 
sired system as specified by the command input file. 

There are two characteristics of the UNAERO 
subroutines that the user must be aware of in order 
to use the subroutines. These characteristics are 
(1) the manner in which the dynamic data storage 
is allocated and (2) the structure of the file for the 
data base containing the description of the vehicle 
that is being investigated. 

The data storage allocation method is similar to 
the method described by Chung-Phillips and Rosen 
in reference 4. Before each subroutine is called, the 
dynamic data storage needed for that subroutine is 
calculated and checked against the amount of storage 
available, which is the length of the data stack. In 
figure 2, the dynamic data storage is obtained from 
the array named STACK and the length of the data 
stack is in the variable named SL. If there is enough 
space in the array STACK for the required storage, 
it is partitioned off for the subroutine and the length 
of the stack is adjusted accordingly and passed to 
the subroutine being called. If the subroutine being 
called also calls another subroutine that requires dy- 
namic data storage, an identical procedure is followed 
for the subsequently called subroutine. On comple- 
tion of the subroutine, the length of the stack is ad- 
justed so the data storage can be used by other sub- 
routines. Since UNAERO has a large number of op- 
tions that affect the amount of dynamic data storage 
required, the storage requirements have never been 
explicitly calculated. The approach that has been 
adopted is to set the length of STACK to a value 
based on experience and run the program. If the 
problem being analyzed requires more than the avail- 
able dynamic storage, the execution of the problem is 
terminated with the appropriate error message. The 
user then adjusts the length of STACK and reruns 
the program. Using the dynamic storage approach 
allows a large aeroelastic problem to be analyzed 
with a modest amount of space. For example, all 
the examples that are presented in this report were 
run with a stack length of 4000 memory locations. 

The other characteristic of UNAERO the user 
must be aware of is the manner in which UN- 
AERO handles the discrete and continuous data 
bases. These files are both stored in the file that 
has its unit designator stored in the variable JD. If 
the user is generating the continuous-frequency data 
base, UNAERO will destroy the discrete-frequency 
data base and replace it with the continuous data 


base being generated. It is the responsibility of 
the user to make sure that the continuous-frequency 
data base is not destroyed when the user saves the 
discrete-frequency data base. On the CDC Network 
Operating System (NOS) this is accomplished by the 
following; 

1. Use an indirect access file to store the continuous 
and discrete data bases. 

2. After the program has executed and the user de- 
sires to retain the continuous-frequency data base, 
store the indirect access file containing the data 

• base under a different name on the permanent 
storage facility. 

The example in figure 2 represents the simplest 
calling sequence of the UNAERO subroutine that 
will calculate a model of the unsteady aerodynamic 
forces. UNAERO requires five data files that are des- 
ignated on the PROGRAM card as INPUT, OUT- 
PUT, VDATA, SCRl, and SCR2. The command 
input file is assigned to the INPUT file, the system 
output file is assigned to the OUTPUT file, and the 
data bases are assigned to the VDATA file. The files 
SCRl and SCR2 are scratch files for temporary stor- 
age. The unit designators for these files are passed to 
UNAERO by the variables IR, IW, JD, JS, and JSl 
in the labeled common block INOUT. The maximum 
dimensions of the system matrices with the vari- 
able names A, B, C, and D are passed to UNAERO 
in the variables NXM, NRM, and NUM. The print 
level and the dynamic storage size are passed to UN- 
AERO with the variables IPRINT and SL. On com- 
pletion of the calculations, UNAERO returns the co- 
efficients of the model of the aeroelastic vehicle, as 
defined by equation (25), in the matrices A, B, C, 
and D. UNAERO returns the actual dimensions of 
the system quadruple as the number of state equa- 
tions defining the model in NX, as the number of 
outputs in NR, and as the number of inputs in NU. 

In the remainder of this section, the calling se- 
quence and parameters of the subroutine UNAERO 
and the variables in the labeled common block IN- 
OUT are listed and their key characteristics are 
discussed. 

Calling sequence: 

CALL UNAERO(A,B,C,D,NXM,NRM,NUM,WKl, 
STACK,SL,NX,NR,NU) 

Argument list: 

A System-matrix dimensioned (NXM x 

NXM) output (see eq. (25)) 

B Control-input-matrix dimensioned 

(NXM X NUM) output (see eq. (25)) 
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C Measurement-matrix dimensioned 

(NRM X NXM) output (see eq. (25)) 

D Control-response-matrix dimensioned 

(NRM X NUM) output (see eq. (25)) 

NXM Maximum number of states in the 

system quadruple (input); number of 
rows and columns of A as dimensioned 
in the calling program 

NRM Maximum number of outputs in the 

system quadruple (input); number of 
rows of C and D as dimensioned in the 
calling program 

NUM Maximum number of inputs in the 

system quadruple (input); number of 
columns of B and D as dimensioned in 
the calling program 

WKl Work array of length NXM 

STACK Work array of length SL (see discus- 
sion on dynamic storage) 

SL Length of array STACK (input) (see 

discussion on dynamic storage) 

NX Actual number of states in the system 

quadruple (output) 

NR Actual number of outputs in the 

system quadruple (output) 

NU Actual number of inputs in the system 

quadruple (output) 

Labeled common block INCUT: 

COMMON/INOUT IR,IW,JD,JS,IPRINT,JS1 

Variable list: 

IR Unit designator for command input 

file 

IW Unit designator for print file 

JD Unit designator for data base file 

JS Unit designator for scratch file 

IPRINT Control variable for print file (see 
table 2 for options) 

JSl Unit designator for a second scratch 

file 

8.2 Description of Subroutines and Functions 

In this section a description of each subroutine is 
given. This section is intended primarily for those 
who are installing and maintaining UNAERO. The 
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information in this section can be skipped if the user 
is interested only in executing the program. The 
structure and hierarchy of the subroutines used in 
UNAERO are given in figure 8. A horizontal line 
between two subroutine names indicates that the 
subroutines on the left call the subroutines on the 
right. In the sections that follow, each of the sub- 
routines and functions in figure 8 is listed in alpha- 
betical order with its usage, purpose, and a descrip- 
tion of its argument list. Exceptions are subroutines 
PRESEN, SENSOR, USER, and CONMIN. Descrip- 
tions of PRESEN and SENSOR can be found in ap- 
pendix A, and an example of subroutine USER can 
be found in appendix B. The subroutine CONMIN 
is described in detail in reference 5. The subrou- 
tine SNVDEC is not part of this development and is 
obtained from reference 6, but it is included in this 
section for completeness. 

8.2.1 Description of global variables. To reduce 
the space used for listing the descriptions of the 
argument lists, the following variable list contains 
the description of the variables that are common to 
several subroutines. The entries for these variables 
have been omitted from the individual subroutine 
argument lists. A variable can be assmned to be an 
input variable unless otherwise noted. 

APLUS matrix of equation (22) and the 

solution matrix of equation (16) 
(dimensioned NPPT x NM). APLUS 
is an output matrix for subroutines 
CALAP and SNVDEC and a work 
matrix for subroutines LSTSOL, 
PADSOL, and PRTDIF. 

CBAR Reference length used to convert 

frequency to nondimensional form 

IFS Array of length NRF defining the 

index to the array RNF. IFS is a work 
array of subroutine CALPR. 

NF Number of structural modes in the 

model 

NM = 2(NRF) 

NP Order of the approximation (i.e., the 

degree of the denominator polyno- 
mial of the unsteady-aerodynamics 
approximation) 

NPPT Order of the numerator polynomial in 
the approximation 

NRF Number of reduced frequency values 


P Vector of length NPPT containing 

the coefficients of the numerator 
polynomial approximation. P is an 
output matrix of subroutines CALPJ 
and LSTSOL. 

Q Complex matrix (dimensioned NF 

X NRF) containing the unsteady- 
aerodynamics data being approxi- 
mated. Q is an output matrix of sub- 
routine PRTDIF. 

QHAT Complex matrix (dimensioned NF x 
NRF) containing the estimate of the 
unsteady aerodynamic forces 

QSNV Work array of length NPPT 

R Vector of length NP containing the 

coefficients of the denominator poly- 
nomial. In the subroutine PADSOL 
this vector contains the initial value 
used for the denominator coefficients 
of the approximation. When subrou- 
tine PADSOL returns to the calling 
program, this vector will contain the 
calculated values of the Pade approxi- 
mation. R is a work array for the sub- 
routine CALPR. 

RNF Array of length NRF containing the 
values of reduced frequency. RNF 
is a scalar containing a value of 
reduced frequency in the functions 
EVALP, EVALR, EVALRR, PHID, 
and PHIVK. 

SIGMA Root-mean-square value of the gust 
(eqs. (23) and (24)) 

SCALE Characteristic length of the gust 
(eqs. (23) and (24)) 

VON, DRY Logical variables specifying which 

scaling routine is called for in the gust 
attenuation option. The routine called 
is defined by the following; 


VON 

DRY 

Routine 

called 

.T. 

.T. 

None 

.T. 

.F. 

PHIVK 

.F. 

.T. 

PHID 

.F. 

.F. 

USER 


VSNV Work array (dimensioned NPPT x 
NPPT) 


WK Work array of length NF 

WKl Work array of length NPPT 

WTF Array of weights (dimensioned NRF) 

Z Complex matrix (dimensioned NF x 

NRF) containing the error between 
the aerodynamics data and approxi- 
mation (see eq. (17)). Z is an output 
matrix of subroutine CALZFQ and a 
work matrix otherwise. 

• 8.2.2 Subroutine CALAP. Subroutine CALAP 
calculates the solution to the least-squares optimiza- 
tion problem described by equation (16). The sub- 
routine forms the matrix A using the selected reduced 
frequencies and calls subroutine SNVDEC to calcu- 
late the pseudoinverse of A. 

Calling sequence: 

CALL CALAP(A,NM,NPPT,NRF,WTF,IFS,RNF, 
APLUS,QSNV,VSNV) 

Argument list: 

A Work array (dimensioned NM x NPPT) 

8.2.3 Subroutine CALGRA. Subroutine CALGRA 
calculates the gradient of the objective function with 
respect to the coefficients of the denominator polyno- 
mial. The Rh element of the vector G is the gradient 
of the ith coefficient of the R polynomial in table 1 
and is calculated with equations (20) to (22). 

Calling sequence: 

CALL CALGRA(G,Z,QHAT,APLUS,WK,WK1,R, 
IFS,RNF,Q,WTF,NP,NPPT, 
NRF,NM,NF) 

Argument list: 

G Output vector of length NP containing 
the gradient of the objective function 

8.2.4 Subroutine CALPJ. This subroutine cal- 
culates the numerator coefficients of the approxima- 
tion by using a least-squares method. This is ac- 
complished by multiplying the pseudoinverse of the 
matrix A (see eq. (16)) by the appropriate vector 
to obtain the solution of the minimization of equa- 
tion (16). The solution is returned in the vector P. 

Calling sequence: 

CALL CALPJ(Q,WTF,WK,R,RNF,IFS,APLUS,P,J, 
NP,NF,NRF,NPPT,NM) 
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Argument list: 

J Integer indicating the structural mode 
number for which the approximation is 
being calculated 

8.2.5 Subroutine CALPR. The subroutine CALPR 
calculates the parameters of the unsteady-aero- 
dynamics approximations for a column of the gener- 
alized aerodynamic-force matrix and selected values 
of the reduced frequency. These parameters are cal- 
culated with the least-squares method or the Fade 
method and are written to the VDATA file. 

Calling sequence: 

CALL CALPR(MS,CS,GS,NF,NC,NG,NFM,NCM, 
NGM,DUM,M1M,M2M,RNF, 
NRFM,IRNF,BUFF,BUF, 
R,WTF,IFS,Q,X,XMAX,CBAR) 


Argument list: 

MS 

Structural mode selection array of 
length NF 

CS 

Control mode selection array of length 
NC 

GS 

Gust mode selection array of length 
NG 

NC 

Number of control modes in the model 

NG 

Number of gust modes in the model 

NFM 

Number of structural modes in the 
discrete data base 

NCM 

Number of control modes in the 
discrete data base 

NGM 

Number of gust modes in the discrete 
data base 

DUM 

Work array (dimensioned MIM x M2M) 

MIM 

= NFM -h NCM 

M2M 

= NFM -t- NCM -t- NGM 

NRFM 

Number of reduced-frequency values in 
the discrete data base 

IRNF 

Work array of length NRFM 

BUFF 

Input data buffer of length 8 

BUF 

Input data buffer of length 20 

X 

Stack array 

XMAX 

Length of stack array 


8.2.6 Subroutine CALZFQ. This subroutine cal- 
culates the estimate of the unsteady-aerodynamic- 
force approximation function (Q in eq. (18)), the er- 
ror vector between the data and the approximation 
{zm in eq. (18)), and the value of the objective func- 
tion that is being minimized {Jj in eq. (19)). 

Calling sequence: 

CALL CALZFQ(QHAT,Q,R,APLUS,WK,WK1, 
NF,NRF,NP,NPPT,NM,Z,F, 
IFS,RNF,WTF) 

Argument list: 

F Output variable F containing the value 
of the objective function 

8.2.7 Subroutine ERMON. This subroutine is 
called when an error in the program is detected. 
The subroutine prints the appropriate error message 
and returns if there is a valid way of continuing and 
generating useful data. Otherwise the subroutine 
terminates the program execution with a CDC NOS 
return code of 105. 

Calling sequence: 

CALL ERMON(ERRNO,PAGE,BUFF) 

Argument list: 

ERRNO Error number being printed (input) 

PAGE Logical input variable that, if true, 

will print the error message at the top 
of the next page and terminate the 
program 

BUFF Input string of up to 80 characters 

used to pass a message to the subrou- 
tine for printing along with the error 
message 

8.2.8 Subroutine ERR. Subroutine ERR is called 
by the CDC CYBER NOS monitor whenever an 
illegal-data-in-field or a data-overflow error is de- 
tected during execution. If the error number passed 
to ERR in the labeled common block ERROR is not 
zero, that error number is passed to the subroutine 
ERMON as the argument ERRNO, which generates 
the appropriate error message. This subroutine is 
used to trap illegal syntax usage in the command. 
For the cases where the command file input is not 
being decoded, the error flag is set to zero so that 
normal detection and processing of data overflow can 
occur. 

8.2.9 Subroutine ERRMON. The subroutine 
ERRMON handles the special error conditions that 
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arise when the maximum system dimensions NXM, 
NRM, and NUM (i.e., the dimensions of the matri- 
ces with the variable names A, B, C, and D in the 
argument list of the subroutine UNAERO) are ex- 
ceeded. When this subroutine is called, the requested 
system dimensions and the maximum system dimen- 
sions are printed and the program is terminated with 
a CYBER NOS return code of 101. 

Calling sequence: 

CALL ERRMON(NX,NC,NR,NXM,NCM,NRM) 
Argument list: 

NX Number of system states being 
requested 

NC Number of control states being 
requested 

NR Number of output states being 
requested 

NXM Maximum number of system states 
NCM Maximum number of control states 
NRM Maximum number of output states 

8.2.10 Function EVALP. The function EVALP 
evaluates the numerator polynomial of the unsteady- 
aerodynamics approximation at the reduced fre- 
quency RNF. 

Calling sequence: 

X=EVALP(P,RNF,NPPT) 

8.2.11 Function EVALR. The function EVALR 
evaluates the denominator polynomial of the 
unsteady-aerodynamics approximation at the re- 
duced frequency RNF. 

Calling sequence: 

X=EVALR(R,RNF,NP) 

8.2.12 Function EVALRR. The function EVALRR 
calculates the value of the partial derivative of the 
denominator polynomial R with respect to the fcth 
coefficient of the polynomial r (see eq. (21)) evalu- 
ated at a reduced-frequency value of RNF. 

Calling sequence: 

X=EVALRR(R,K,RNF,NP) 


Argument list: 

K Integer indicating the parameter that 
the partial derivative is to be taken to 

8.2.13 Subroutine FABCD. Subroutine FABCD 
forms the system quadruple, matrices A, B, C, and 
D, from the data stored in the continuous-frequency 
data base. 

Calling sequence: 

CALL FABCD(A,B,C,D,XN1,XN2,NXM,NRM, 
NUM,CBAR,V,QBAR,NF,NC, 
NG,NS,NGF,XN3,XN4,XN5, 
NGFA,XN6,XN7,XN8,XN9, 
XN10,XN11,XN12,XN13, 
XN14,NX,NR,NU) 

Argument list: 

A System matrix defined by equa- 

tion (25) (dimensioned NXM x 
NXM) 

B System matrix defined by equa- 

tion (25) (dimensioned NXM x 
NUM) 

C System matrix defined by equa- 

tion (25) (dimensioned NRM x 
NXM) 

D System matrix defined by equa- 

tion (25) (dimensioned NRM x 
NUM) 

XN1,XN9, Work arrays (dimensioned 

XNll NXM X NXM) 

XN2 Work array (dimensioned 

NRM X NXM) 

NXM Maximum number of states that the 

system matrices can accommodate, 
as defined by the calling program 

NRM Maximum number of outputs 

that the system matrices can 
accommodate, as defined by the 
calling program 

NUM Maximum number of inputs that 

the system matrices can accom- 
modate, as defined by the calling 
program 

V Reference velocity 

QBAR Reference dynamic pressure 

NC Number of control modes in the 

model 
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NG 

Number of gust modes in the model 

NS 

Number of system outputs 

NGF 

Number of generalized forces to be 
retained 

XN3 

Work array of length 4 

XN4 

Work array (dimensioned NF x 7) 

XN5 

Work array (dimensioned NF x 8) 

NGFA 

Vector of length NGF containing an 
index as to which generalized forces 
are to be retained in the system 
being formed 

XN6,XN7, 

XN8 

Work arrays (dimensioned 
NF X NF) 

XNIO 

Work array of length NXM 

XN12,XN13, 

XN14 

Work array (dimensioned NS x NF) 

NX 

Number of system states generated 
by this subroutine 

NR 

Number of system outputs gener- 
ated by this subroutine 

NU 

Number of system inputs generated 
by this subroutine 


8.2.14 Subroutine FDATA. The subroutine 
FDATA reads and modifies the data from the 
discrete-frequency data base and returns the data in 
the complex matrix Q of equation (18) according to 
the mode selection array and frequency selection ar- 
ray in the argument list. Before returning to the 
calling program the data are attenuated if the user 
has so specified, specified. 

Calling sequence: 


NRFM Number of reduced-frequency values in 
the discrete-frequency data base 

IFORMl Format for reading data from the 

discrete-frequency data base (VDATA) 

8.2.15 Subroutine LSTSOL. Subroutine LSTSOL 
calculates the least-squares solution for the numer- 
ator coefficients of the unsteady-aerodynamics ap- 
proximation and writes the results in the continuous- 
frequency data base. 

Calling sequence: 

CALL LSTSOL(A,B,NM,NRF,Q,NF,WTF,IFS, 
RNF,R,NP,P,NPPT,QSNV, 
VSNV,APLUS,IFORM) 

Argument list: 

A Work array (dimensioned NM x 

NPPT) 

B Work array of length NM 

IFORM Format for writing the approximation 
to the continuous-frequency data base 

8.2.16 Subroutine PADSOL. The subroutine 
PADSOL calculates the Pad6 coefficients of the 
unsteady-aerodynamics approximation and writes 
the results in the continuous-frequency data base. 

Calling sequence: 

CALL PADSOL(R,G,VLB,VUB,SCAL,DF,AX,S, 
G1,G2,B,C,ISC,IC,MS1X,N1P, 
N2P,N3P,N4P,N5P,NP,Q, 
APLUS,Z,QHAT,WK,WK1, 
IFS,RNF,WTF,NPPT,NRF,NM, 
NF,IFORM,VSNV,A,AT) 

Argument list: 


CALL FDATA(DUM,M1M,M2M,Q,NF,NRF,IFS, 

G,VLB, 

Work arrays of length N IP 


MS, IM, NRFM, IFORM 1 ,SIGM A, 

VUB,SCAL, 

used by the subroutine CONMIN 


SCALE,RNF,CBAR,VON,DRY) 

DF,S 


Argument list: 

AX 

Work array (dimensioned NIP 

DUM 

Work array (dimensioned MIM x M2M) 


X N3P) used by the subroutine 
CONMIN 

MIM 

= NF -1- NC, input as defined in 
subroutine CALPR 

G1,G2,ISC 

Work arrays of length N2P used by 
the subroutine CONMIN 

M2M 

= NF -f NC -1- NG, input as defined 
in subroutine CALPR 

B 

Work array (dimensioned N3P 
X N3P) used by the subroutine 

MS 

Mode selection array 


CONMIN 

IM 

Mode for which data are being 

C 

Work array of length N4P used by 


requested 


the subroutine CONMIN 
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IC Work array of length N3P used by 

the subroutine CONMIN 

MS IX Work array of length N5P used by 

the subroutine CONMIN 

NIP = NP + 2 

N2P = 2(NP) 

N3P = NP + 1 

N4P = NP + 1 

N5P = 2(NP + 1) 

IFORM Format for writing the calculated 

approximation to the continuous- 
frequency data base 

A Work array (dimensioned 

NM X NPPT) 

AT Work array (dimensioned NF x 

NPPT); may be made equivalent to 
the matrix A in the calling program 

8.2.17 Function PHID. The function PHID re- 
turns the attenuation factor for a Dryden turbulence 
spectrum (eq. (24)) at the spectral frequency that is 
equivalent to the reduced frequency specified by the 
parameter RNF. 

Calling sequence: 

X=PHID(SIGMA, SCALE, RNF, CBAR) 

8.2.18 Function PHIVK. The function PHIVK 
returns the attenuation factor for a Von Karman tur- 
bulence spectrum (eq. (23)) at the spectral frequency 
that is equivalent to the reduced frequency specified 
by the parameter RNF. 

Calling sequence: 


Argument list: 

AR Input matrix to be printed (dimen- 
sioned ND X NC) 

ND First dimension of NR in the calling 
program 

NR Number of rows to be printed 
NC Number of columns to be printed 

8.2.20 Subroutine PRTDIF. The subroutine 
PRTDIF prints program output as described in 
section 7. 

Calling sequence: 

CALL PRTDIF(QHAT,Q,R,APLUS,WK,WK1, 
NF,NRF,NP,NPPT,NM,Z,F,IFS, 
RNF,WTF,IMX,IM) 

Argument list: 

F On return, contains the value of 

the objective function for the given 
parameters 

IMX Mode number in the discrete- 
frequency data base 

IM Mode number in the continuous- 
frequency data base 

8.2.21 Subroutine REFORM. The subroutine 
REFORM converts the coefficients of an unsteady- 
aerodynamics approximation from the form that is 
used in the calculations (eqs. (3) and (4)) into the 
form used by the subroutine FABCD (eq. (5)). 

Calling sequence: 

CALL REFORM(P,R,A2,A1,AO,DO,D1,D2,D3,W, 
NP,NF,CBAR,V) 


X=PHIVK(SIGMA,SCALE,RNF,CBAR) 

8.2.19 Subroutine PRNT. The subroutine PRNT 
is a matrix utility that prints large matrices on 
the output device designated by the variable IW in 
the labeled common block INOUT. The subroutine 
divides the large matrix into submatrices that will fit 
on an 8.5 x 17-in. page. Each of the submatrices is 
then printed. 

Calling sequence: 

CALL PRNT(AR,ND,NR,NC) 


Argument list: 

A2,A1,A0, Work array of length NF 
DO, Dl, 

D2,D3 

W Work array of length NP 

V Reference velocity 

8.2.22 Subroutine SNVDEC. The subroutine 
SNVDEC calculates the pseudoinverse of a matrix. 
Additional information can be found in reference 6. 
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Calling sequence: 

CALL SNVDEC(IOP,MD,ND,M,N,A,NOS,B,IAC, 
ZTEST,Q,V,IRANK,APLUS,IERR) 

Argument list: 

lOP Set equal to 4 for pseudoinverse option 

MD Input integer specifying the maximum 

row dimension for A 

ND Input integer specifying the maximum 

row dimension for V 

M Input integer specifying the number of 

rows in A 

N Input integer specifying the number of 

columns in A 

A Input matrix, destroyed by the 

program 

NOS Number of columns in B. The value of 
zero is used by the calling program. 

B Dummy argument 

lAC Input integer specifying the number 

of decimal digits of accuracy in the 
elements of the matrix A 

ZTEST Output variable used by SNVDEC to 
test for zero 

Q Vector of length N which will contain 

the singular values on return from the 
subroutine 

V Work array (dimensioned ND x N) 

IRANK Rank of the matrix A 

lERR Returned error indicator, as follows; 

= 0, no error 

= K, Kth singular value not found 
after 30 iterations 

= — 1, this matrix is close to a matrix 
which is of lower rank than IRANK 

8.2.23 Subroutine ZZLINE. The subroutine 
ZZLINE formats the printed output into pages of 
44 lines each. If the number of lines N being added to 
the output file is greater than 44, a page eject com- 
mand is issued and the line counter is set to N. If N is 
greater than 44, then the line counter is set to zero. 
The unit designator for the output is contained in the 
variable IW of the labeled common block INOUT. 


Calling sequence: 

CALL ZZLINE(N) 

Argument list: 

N Number of lines to be added to the 
current page 

8.2.24 Subroutine ZZPAGE. The subroutine 
ZZPAGE is used with the subroutine ZZLINE to re- 
turn 44, the number of lines printed since the last 
page eject. 

Calling sequence: 

CALL ZZPAGE(N) 

Argument list: 

N Number of lines left on the page 

9 Examples of UNAERO Use 

Several examples are presented which illustrate 
the options that are available in the program and 
which show the input used to invoke these options. 
The output listings that are produced by running 
the program with these input are omitted. This 
output is omitted in the interest of saving space and 
has the same form as described in section 7, but it 
has different data printed for the results. All the 
examples presented use the same task I discrete- 
frequency data base as listed in figure 5. This data 
base is described in the next paragraph followed by 
a discussion of all the examples. 

The discrete-frequency data bases for the exam- 
ples described are based on the aeroelastic semispan 
wind tunnel model (ARW-1) that was used for the ex- 
ample in reference 1. The wing geometry is presented 
in figure 9 and the generalized masses and frequencies 
which describe the aeroelastic model are presented in 
table 3. Structural damping is assumed to be negli- 
gible. The simple harmonic oscillatory aerodynamic 
forces are calculated with a doublet-lattice technique 
similar to that described in reference 7. In order to 
calculate the pressure distribution on an oscillating 
wing undergoing simple harmonic motion, the lift- 
ing surface is subdivided into an array of trapezoidal 
boxes arranged in strips parallel to the air stream, as 
shown in figure 10. The lifting surface is then rep- 
resented by a lattice of doublets located along the 
quarter-chord of each box. The downwash condi- 
tion is satisfied at the three-quarter-chord of each 
box by equating it to the downwash resulting from 
the slope and deflection rate of each structural mode. 
The semispan lifting surface is divided into 210 boxes 
arranged in 30 spanwise strips of 7 chordwise boxes. 
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The side-edge effects of the control surface (as de- 
scribed in ref. 7) are not applied. Oscillatory aero- 
dynamic forces are calculated at eight reduced fre- 
quencies (0, 0.1, 0.3, 0.5, 0.7, 0.9, 1.3, and 1.8). The 
structural model consists of 10 structural modes, 1 
control-surface mode, and 1 gust mode. 

9.1 Example 1 

The command input file for example 1 is pre- 
sented in figure 1(a). In this example a first-order 
Fade approximation is used for 5 of the 10 structural 
modes that are stored in the discrete data base. The 
modes that are retained are 1, 2, 4, 5, and 6. (For 
the ARW-1 wind tunnel model, the third mode is 
mainly motion fore and aft in the plane of the wing 
and contributes little to the dynamic behavior of the 
wing under investigation.) Since the gust input mode 
is not needed for the flutter calculation, it is also 
deleted. This illustrates how to delete a set of modes 
from the data being used for the calculation of the 
approximation. Each approximation is based on data 
calculated at the first six values of reduced frequency. 
The last six lines of figure 1 (a) represent a second call 
to subroutine UNAERO. In this call, UNAERO uses 
the stored data in the continuous-frequency data base 
and calculates a system quadruple (i.e., the matrices 
A, B, C, and D) at a new set point with the dynamic 
pressure increased by a factor of two. 

9.2 Example 2 

The command input file for example 2 is pre- 
sented in flgure 1(b). This example is identical to the 
first example except that now second-order Fade ap- 
proximations are being calculated. The initial value 
of the denominator coefficients is based on the results 
from example 1. One initial value is set to a fraction 
of the result generated from example 1 and the other 
initial value is set to a multiple of that result. For 
this example the initial values are 1/2 and 2 times 
the result of example 1. 

9.3 Example 3 

Example 3, presented in figure 1(c), is also a 
derivative of example 1, except that here fourth-order 
least-squares approximations are used. The least- 
squares approximations are used in this case after 
the denominator coefficients are determined with the 
Fade approximation. 

Using the least-squares approximations with the 
results from the Fade approximations effects a sav- 
ings in computer resources while giving identical re- 
sults. This is often done while parametric studies 
are performed that do not influence the aerodynam- 
ics approximation calculations. 


9.4 Example 4 

Figure 1(d) is the command input file which il- 
lustrates the generation of a system quadruple with 
a weighted least-squares approximation. The system 
generated is of “minimum order” in that the zero- 
order least-squares approximations induce no addi- 
tional states into the system to model the effects 
of the unsteady aerodynamics. This approximation 
function results in the largest errors of any of the 
approximation functions that UNAERO can gener- 
ate. With proper selection of the weights, however, 
useful results can be obtained. In fact, this exam- 
ple predicts the flutter dynamic pressure to the same 
degree of accuracy as does the second-order Fade ap- 
proximation, since a priori knowledge was used in the 
selection of the weight so that the minimum error oc- 
curs at the reduced frequency of flutter. The weights 
are selected to be inversely proportional to the dif- 
ference between the reduced frequency of the data 
and the reduced frequency of flutter. Note that the 
weights for the data at the seventh and eighth values 
of reduced frequency have been left blank, thereby 
resulting in a zero for the weight at those frequen- 
cies. Using a weight of zero requires more computer 
resources than does removing that value of reduced 
frequency from the data base before the approxima- 
tion is calculated. In this example, the frequencies 
were retained, since it was desired to observe the cal- 
culated value of the approximation at these reduced 
frequencies. 

9.5 Example 5 

The conunand input file for this example is pre- 
sented in figure 1(e). In this example, two struc- 
tural modes (modes 1 and 2) and the gust input 
are retained whereas the control-surface mode is re- 
moved. The approximations to the unsteady aero- 
dynamic forces due to motion in the first mode 
are fourth-order Fade approximation functions based 
on the first six frequencies in the discrete data 
base. The unsteady-aerodynamic-force approxima- 
tions corresponding to the second mode are fourth- 
order least-squares approximations, but these ap- 
proximations are based on all the frequencies in the 
discrete data base. The gust input is approximated 
by Fade approximations with second-order numera- 
tors and third-order denominators. These approxi- 
mations are based on data at three values of reduced 
frequency (0, 0.3, and 0.9) and are attenuated with 
a Von Karman gust spectrum with a characteristic 
length of 1200 in. and a vertical velocity of 12 in/sec. 
Two sensors are to be added to the system equations 
with the user-supplied FRESEN and SENSOR sub- 
routines. An explanation of these two subroutines is 
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given in appendix A. The subroutine SENSOR reads 
the data for the location of the sensors from the two 
cards after the “ADD 00002 SENSORS AT” card and 
provides UNAERO with the required information to 
calculate the sensor equations. The reader may refer 
to appendix A for any further discussion. 

10 Concluding Remarks 

UNAERO is a collection of FORTRAN IV sub- 
routines for approximating the unsteady aerody- 
namic forces in the time domain and obtaining a set 
of constant-coefficient first-order differential equa- 
tions. This report serves as an instruction and main- 
tenance manual for UNAERO. These subroutines 
are intended to be run on CDC CYBER computer 
systems. 

UNAERO has provisions for adjusting the num- 
ber of modes used for calculating the approximations 


so that an accurate approximation is generated. The 
number of data points at different values of reduced 
frequency can also be varied to adjust the accuracy of 
the approximation over the reduced-frequency range. 
The denominator coefficients of the approximation 
may be calculated by means of a gradient method or 
may be entered into a least-squares approximation 
technique. Both approximation methods use weights 
on the residual error. UNAERO uses data bases to 
store the coefficients of the approximation. In do- 
ing so, a new set of system equations, at a different 
dynamic pressure, can be generated without the ap- 
proximations being recalculated. 

NASA Lsingley Research Center 
Hampton, VA 23666-5225 
April 29, 1985 



Appendix A 

Sample Subroutines PRESEN and 
SENSOR 

The mode shapes that are used to describe the 
dynamics of the vehicle are normally provided as a 
set of modal displacements at specific points on the 
vehicle. The generation of the sensor equations that 
define the output of a sensor (at an arbitrary location 
on the surface of the vehicle) requires that the value 
of the displacement for each mode be calculated at 
the location of the sensor. This calculation requires 
the use of an interpolation method for the modes 
and can depend on the structural model used in 
generating the modes. The inclusion of all known 
useful interpolation schemes was beyond the scope 
of this work. In order to have the ability to add 
sensor equations to the system equations and to meet 
the individual needs of the users, the section of the 
program which performs the modal interpolation is 
left for the user to provide. 

UN AERO assumes that the output of a sensor 
can be modeled as a linear combination of the ac- 
celeration, the velocity, and the displacement of the 
generalized coordinates. Thus the sensor output is 
modeled by the following equation: 

Y = AP * q -t- VP * q -t- DP * q (Al) 

where Y is a vector of the sensor outputs, q is a 
vector of the generalized coordinates, and AP, VP, 
and DP are matrices of the coefficients which relate 
the generalized coordinate’s acceleration, velocity, 
and displacement to the sensor output. Dots over 
q denote a derivative with respect to time. 

To interface the interpolation method UNAERO 
requires that the user provide two subroutines. One 
subroutine, PRESEN, is called during the execution 
of task I. This subroutine is required to edit and 
copy the structural modal data, which have been 
appended by the user to the end of the discrete- 
frequency data base, in the continuous-frequency 
data base. The other subroutine, SENSOR, is called 
during task II. It uses the copied data from the 
continuous-frequency data base to calculate the re- 
quired coefficients of the sensor equation. 

In order for the subroutine PRESEN to edit the 
discrete-frequency data base, an array of integers 
indicating the modes which are to be retained for 
the approximations is transferred to the subroutine 
in the variable array MS. Each element of MS is an 
integer that is a mode number to be retained. The 
numbers in MS are initialized by the “mode selection 
array” in figure 4(d). If the user does not reduce 


the number of modes used for the approximation, 
the MS array will contain the numbers 1 to NFM. 
If the user does omit some of the modes, the array 
will be of length NF and contain the “mode selection 
array” data that was read by the program from the 
command file. 

The variables that contain the INPUT/OUTPUT 
unit designators for the continuous- and discrete- 
frequency data bases are in the labeled common block 
INOUT. The INPUT/OUTPUT unit designator for 
the discrete-frequency data base is contained in the 
variable JSl, and the INPUT/OUTPUT unit desig- 
nator for the continuous-frequency data base is con- 
tained in the variable JD. On entry into the sub- 
routine PRESEN, the discrete data base must be re- 
wound and positioned at the beginning of the modal 
displacement data. 

The calling sequence for the subroutine is 
CALL PRESEN(NFM,NF,MS,IFORM,IFORMl) 
where 

NFM Number of structural modes in the 
discrete-frequency data base 

NF Number of structural modes in the 

continuous-frequency data base 

MS Array of length NF. Each element of 

MS contains the number of the mode 
that is to be copied from the discrete 
data base to the continuous data base. 

IFORM Format for writing data to the contin- 
uous data base 

IFORM 1 Format for reading data from the 
discrete data base 


A sample of the subroutine PRESEN is presented in 
figure Al. 

Subroutine SENSOR is called during task II by 
UNAERO to calculate the elements of the AP, VP, 
and DP matrices in equation (Al). Each of these ma- 
trices has dimensions of NS x NF, where NS is the 
number of sensors (as specified in the “ADD SEN- 
SORS” command) and NF is the number of modes. 
For example, if a single accelerometer is added to 
a system with four modes, the matrices in equa- 
tion (Al) would have dimensions of one by four in 
subroutine SENSOR. The subroutine sets each ele- 
ment of the AP matrix to the interpolated modal dis- 
placement for each mode. On entry to the subroutine 
SENSOR, the continuous-frequency data base is po- 
sitioned to the first card of the augmented data that 
was written to the data base by the subroutine PRE- 
SEN. The INPUT/OUTPUT unit designator for the 
continuous-frequency data base is contained in the 
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variable JS in the labeled common block INOUT. 
The command input file, with the input/output unit 
designator contained in the variable JD of the la- 
beled common block INOUT, is positioned to the 
next record after the “ADD SENSOR AT” request. 
The user can use input following this card to pass 
data and parameters, such as the sensor position, to 
the interpolation routine contained in the subroutine 
SENSOR. 


The calling sequence for the subroutine SENSOR 
is 

CALL SENSOR(AP,VP,DP,NS,NF,IFORM) 

where 


AP 

Output matrix of coefficients that 
relates the output of the sensor to 
the acceleration of the generalized 
coordinate 

VP 

Output matrix of coefficients that 
relates the output of the sensor to the 
velocity of the generalized coordinate 

DP 

Output matrix of coefficients that 
relates the output of the sensor to the 
position of the generalized coordinate 

NS 

Number of sensors being added 

NF 

Number of generalized coordinates in 
the system 

IFORM 

Format for reading the discrete data 
base 


The subroutines PRESEN and SENSOR are re- 
quired only when the user requests that sensors be 
added to the equations being generated. If the user 
makes the request “ADD SENSORS AT” in the 


command input file without supplying the subrou- 
tines PRESEN and SENSOR, an error message is 
generated. 

C SAMPLE PRESEN SUBROUTINE 

SUBROUTINF, PRESEN( NFM ,NF ,MS , IFORM , IFORMl ) 
DIMENSION MS(NF),BUF(8),Z(102,10) 

COMMON /INOUT/ IR,IW,JD,JS,IPRINT,JS1 
C REWIND DISCRETE DATA FILE 

REWIND JSl 
HNODE=10HNODES 

1 FORMAT(BAIO) 

2 F0RMAT(I5) 

3 FORMAT(25HZ - DISPLACEMENT FOR MODE ,15) 

C SEARCH FOR BEGINNING OF MODEL DATA 

100 CONTINUE 

READ(JS1,1) TEMP 
IF(TEMP.EQ.HNODE) GO TO 110 
GO TO 100 
no CONTINUE 

C READ NUMBER OF NODES 

READ(JS1 ,2) NN 

C READ AND WRITE COMMENT CARD FOR X-NODES 

READ(JS1,1) BUF 
WRITE(JD,1) BUF 

C READ AND WRITE X-NODES 

READ( JSl, IFORMl) (Z(I,1) ,1-1 ,102) 

WRITE! JD, IFORM) (Z( I , 1 ) , I-l , 102) 

C READ AND WRITE COMMENT CARD FOR Y-NODES 

read! JS l, 1) BUF 
WRITE!JD,1) BUF 

C READ AND WRITE Y-NODES 

READ! JSl , IFORMl) !Z!l , 1 ) , 1=1 , 102 ) 

WRITE! JD, IFORM) !Z!l , 1 ) , 1=1 , 102) 

C LOOP FOR EACH MODE 

DO 200 1=1,10 

C SKIP HEADER AND READ Z-DISPLACEMENT FOR MODE I 

READ! JSl) 

READ! JSl , IFORMl) !Z! J ,1) , J=1 , 102) 

200 CONTINUE 

C LOOP FOR EACH OUTPUT MODE 

DO 210 1=1 ,NF 

C WRITE HEADER 

WRITE!JD,3) I 

C WRITE Z-DISPLACEMENTS FOR SELECTED MODE 

WRITE! JD, IFORM) !Z!J,MS!D) ,J=1 ,102) 

210 CONTINUE 
RETURN 
END 

Figure Al. Sample of PRESEN subroutine. 
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Appendix B 
Subroutine User 

It is recognized that the user may wish to ap- 
ply modifications to the data before the unsteady- 
aerodynamics approximations are calculated. This 
subroutine allows the user to apply a change in the 
magnitude of the data as a function of the frequency. 
A possible use for this change is to attenuate the data 
for the gust mode with a gust spectrum which is nei- 
ther a Von Karman nor a Dryden distribution. This 
might occur in an experimental test with the random 
disturbance being artificially generated. Subroutine 
USER is called instead of the function PHIVK or 
PHID whenever the “ATTENUATE DATA” card 
specifies “USER SPECTRUM.” (See fig. 4(k).) 

The calling sequence for the subroutine USER is 

X=USER(SIGMA, SCALE, RNF, GEAR) 


The parameters SIGMA, SCALE, and GEAR are 
read in from the command input file. For this option 
the parameters SIGMA and SCALE pass to the 
function USER any data the user may desire. The 
parameter RNF is the reduced frequency for which 
the calling program (subroutine FDATA) requests 
the attenuation factor. The parameter GEAR is 
the reference length from the input described in 
figure 4(f). The formula used to attenuate the data 
is 

(NEW DATA) = (OLD DATA)*USER(SIGMA, 
SCALE,RNF,CEAR) 

where (NEW DATA) is the data used to calculate the 
approximation and (OLD DATA) is the correspond- 
ing value in the discrete data base. 
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TABLE 1. PARAMETERS OF DENOMINATOR POLYNOMIAL R{joj) 


Order of 
approximation 

Polynomial used in equation (4) 

1 

(jk) + ri 

2 

+ r2{jk) + ri 

3 

[(jk)^ + r 2 {jk) + ri][(yfc) + r^] 

4 

[{ifc)2 + r2{jk) + ri][(jfc)2 + r^ijk) + ra] 


TABLE 2. PRINT OPTIONS 


IPRINT 

Resulting output 

0 

Print nothing. 

2 

Print command input file (item 1, fig. 7). 

4 

Print command input file, parameter of the 
approximation, original data, approximated 
data, and error and system matrices 
(items 1-13, fig. 7). 

6 

Print diagnostic information plus the 
information of IPRINT=4. 

7 

Set CONMIN to print level 6 plus the 
information of IPRINT=6. (See microfiche 
output.) 


TABLE 3. FREQUENCIES AND GENERALIZED MASSES 
OF ARW-1 WIND TUNNEL MODEL 


Mode 

Natural frequency, 
Hz 

Generalized mass, 
kg 

1 

5.233 

3.678 

2 

19.129 

7.769 

3 

20.906 

7.044 

4 

25.769 

2.970 

5 

46.110 

4.714 

6 

61.234 

4.758 

7 

79.682 

5.156 

8 

86.030 

11.297 

9 

98.087 

7.558 

10 

118.150 

5.501 
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UNSTEADY AERO DATA 
REFERENC E LE NGTH 25.4 
REFERENC E VE LOCITY 5400. 

REFERENC E DYN AMIC PRESSURE 100. 

USE 00005 ST RUCT URAL MODES 
1 2 4 5 6 

00000 GU ST I NPUTS 

USE A 1 ORDER PADE ' APPROXIMATION FOR MODE 1 
1.02366 

BASE THE APPROXIMATION ON 00006 FREQUENCIES 
1 2 3 4 5 6 

USE A 1 ORDER PADE' APPROXIMATION FOR MODE 2 
1.2135 

BASE THE APPROXIMATION ON 00006 FREQUENCIES 
1 2 3 4 5 6 

USE A 1 ORDER PADE' APPROXIMATION FOR MODE 3 
1.1353 

BASE THE APPROXIMATION ON 00006 FREQUENCIES 
1 2 3 4 5 6 

USE A 1 ORDER PADE' APPROXIMATION FOR MODE 4 
1.1281 

BASE THE APPROXIMATION ON 00006 FREQUENCIES 
1 2 3 4 5 6 

USE A 1 ORDE R PA DE' APPROXIMATION FOR MODE 5 
1.09640 

BASE THE APPROXIMATION ON 00006 FREQUENCIES 
1 2 3 4 5 6 

USE A 1 ORDER PADE' FOR CONTROL 1 
1.3336 

BASE THE APPROXIMATION ON 00006 FREQUENCIES 
1 2 3 4 5 6 

END OF UNSTEADY AERO DATA 
UNSTEADY AERO DATA 
USE APPROXIMATION ON VDATA FILE 
REFERENC E LE NGTH 25.4 
REFERENC E VEL OCITY 5400. 

REFERENC E DY NAMIC PRESSURE 200. 

END OF UNSTEADY AERO DATA 


(a) Command input file for example 1. Key elements are underlined. 
Figure 1. Examples of command input files. 



UNSTEADY AERO DATA 
REFERENCE LENGTH 25.4 
REFERENCE VELOCITY 5400. 

REFERENCE DYNAMIC PRESSURE 100. 

USE 00005 STRUCTURAL MODES 
1 2 4 5 6 

00000 GUST INPUTS 

USE A 2 ORDER PADE' APPROXIMATION FOR MODE 1 
.126136 .50226 

BASE THE APPROXIMATION ON 00006 FREQUENCIES 
1 2 3 4 5 6 

USE A 2 ORDER PADE* APPROXIMATION FOR MODE 2 
.04251 .29160 

BASE THE APPROXIMATION ON 00006 FREQUENCIES 
1 2 3 4 5 6 

USE A 2 ORDER PADE* APPROXIMATION FOR MODE 3 
.02042 .20211 

BASE THE APPROXIMATION ON 00006 FREQUENCIES 
1 2 3 4 5 6 

USE A 2 PADE* ORDER APPROXIMATION FOR MODE 4 
.02613 .22861 

BASE THE APPROXIMATION ON 00006 FREQUENCIES 
1 2 3 4 5 6 

USE A 2 ORDER PADE* APPROXIMATION FOR MODE 5 
.03369 .25956 

BASE THE APPROXIMATION ON 00006 FREQUENCIES 
1 2 3 4 5 6 

USE A 2 ORDER PADE* FOR CONTROL 1 
.15572 .55806 

BASE THE APPROXIMATION ON 00006 FREQUENCIES 
1 2 3 4 5 6 

END OF UNSTEADY AERO DATA 


(b) Command input file for example 2. 
Figure 1. Continued. 



UNSTEADY AERO DATA 
REFERENCE LENGTH 25.4 
REFERENCE VELOCITY 5400. 

REFERENCE DYNAMIC PRESSURE 100. 

USE 00005 STRUCTURAL MODES 
1 2 4 5 6 

00000 GUST INPUTS 

USE A 4 ORDER LEAST SQUARES APPROXIMATION FOR MODE 1 
.23263 1.9581 .50516 .68992 

BASE THE APPROXIMATION ON 00006 FREQUENCIES 
1 2 3 4 5 6 

USE A 4 ORDER LEAST SQUARES APPROXIMATION FOR MODE 2 
.34664 3.2804 .15091 .74886 

BASE THE APPROXIMATION ON 00006 FREQUENCIES 
1 2 3 4 5 6 

USE A 4 ORDER LEAST SQUARES APPROXIMATION FOR MODE 3 
1.9397 3.1545 .033269 .40409 

BASE THE APPROXIMATION ON 00006 FREQUENCIES 
1 2 3 4 5 6 

USE A 4 ORDER LEAST SQUARES APPROXIMATION FOR MODE 4 
1.2788 2.2638 .052972 .48044 

BASE THE APPROXIMATION ON 00006 FREQUENCIES 
1 2 3 4 5 6 

USE A 4 ORDER LEAST SQUARES APPROXIMATION FOR MODE 5 
.57787 1.1957 .058564 .49398 

BASE THE APPROXIMATION ON 00006 FREQUENCIES 
1 2 3 4 5 6 

USE A 4 ORDER LEAST SQUARES FOR CONTROL 1 
.30149 3.61 1.853 3.61 

BASE THE APPROXIMATION ON 00006 FREQUENCIES 
1 2 3 4 5 6 

END OF UNSTEADY AERO DATA 


(c) Command input file for example 3. 
Figure 1. Continued. 



UNSTEADY AERO DATA 
REFERENCE LENGTH 25.4 
REFERENCE VELOCITY 5400. 

REFERENCE DYNAMIC PRESSURE 100. 

USE 00004 STRUCTURAL MODES 
12 4 5 

00000 GUST INPUTS 

USE A 0 ORDER LEAST SQUARES APPROXIMATION FOR MODE 1 
WEIGHT EACH FREQUENCY WITH 

69.44 2500. 30.86 6.93 2.97 

USE A 0 ORDER LEAST SQUARES APPROXIMATION FOR MODE 2 
WEIGHT EACH FREQUENCY WITH 

69.44 2500. 30.86 6.93 2.97 

USE A 0 ORDER LEAST SQUARES APPROXIMATION FOR MODE 3 
WEIGHT EACH FREQUENCY WITH 

69.44 2500. 30.86 6.93 2.97 

USE A 0 ORDER LEAST SQUARES APPROXIMATION FOR MODE 4 
WEIGHT EACH FREQUENCY WITH 

69.44 2500. 30.86 6.93 2.97 

USE A 0 ORDER LEAST SQUARES FOR CONTROL 1 
WEIGHT EACH FREQUENCY WITH 

69.44 2500. 30.86 6.93 2.97 

END OF UNSTEADY DATA 


1.64 


1.64 


1.64 


1.64 


1.64 


(d) Command input file for example 4. 
Figure 1. Continued. 


UNSTEADY AERO DATA 
REFERENCE LENGTH 25.4 
REFERENCE VELOCITY 5400. 

REFERENCE DYNAMIC PRESSURE 100. 

USE 00002 STRUCTURAL MODES 
1 2 

00000 AERODYNAMIC CONTROLS 

USE A 4 ORDER FADE' APPROXIMATION FOR MODE 1 
.08 .6 .48 1.6 

BASE THE APPROXIMATION ON 00006 FREQUENCIES 
1 2 3 4 5 6 

USE A 4 ORDER LEAST SQUARES APPROXIMATION FOR MODE 2 
.34664 3.2804 .15091 .74886 

USE A 2,3 ORDER FADE' APPROXIMATION FOR GUST 1 
1.5 7.5 .1 

BASE THE APPROXIMATION ON 00003 FREQUENCIES 
1 3 6 

ATTENUATE DATA WITH VON KARMAN SPECTRUM 
SCALE 1200. 

RMS 12. 

ADD 00002 SENSORS AT 
73.1 65.0 

47.8 43.5 

END 


(e) Command input file for example 5. 
Figure 1. Concluded. 



PROGRAM EXAMPLE( INPUT, OUTPUT, VDATA, SCRl ,SCR2, 

TAPES =I NPUT , TAPE6 =OUTPUT , TAPE 1 0 =VDATA , 
TAPE12=SCR1 ,TAPE13=SCR2) 

INTEGER SL 

DIMENSION A(32,32),B(32,5),C(7,32),D(7,5), 
STACK(4000),WK02) 

COMMON /INOUT/ IR,IW,JD,JS,IPRINT,JS1 

IR=5 

IW=6 

JD=10 

JS=12 

JS1=13 

NXM=32 

NRM=7 

NUM=5 

IPRINT=4 

SL=4000 

CALL UNAERO(A,B,C,D,NXM,NRM,NUM,WK,STACK,SL,NX,NR,NU) 

STOP 

END 


Figure 2. Sample calling program which uses UNAERO. 


1 2 3 4 5 

12345678901234567890123456789012345678901234567890 
REFERENC E LE NGTH 10.0 
REFERENC E VE LOCITY 120. 

REFERENC E DY NAMIC PRESSURE 100. 

USE A PPROXIMATION ON VDATA FILE 
USE 00002 S TRUCT URAL MODES 
00001 AE RODY NAMIC CONTROLS 
00003 GU ST IN PUTS 

USE A 1 ORDE R LEA ST SQUARES APPROXIMATION 

USE A 2 ORDE R PA PE' APPROXIMATION 

USE A 1,2 ORDER LEAS T SQUARES APPROXIMATION 

USE A 2,3 ORDER PAPE ' APPROXIMATION 

BASE THE APPROXIMATION ON 00004 FREQUENCIES 

WEIG HT EACH FREQUENCY WITH 

RETA IN 00002 NORMALIZED GENERALIZED FORCES AS OUTPUTS 
ADD 00003 SENSORS AT 

ATTENUATE DATA WITH VON KARMAN SPECTRUM 
ATTENUATE DATA WITH DRYD EN SPECTRUM 
ATTENUATE DATA WITH USER SPECTRUM 
SCALE 4. 

RMS 100. 

END 


Figure 3. UNAERO legal command inputs. Key elements are underlined. 






























(e) Statements for using previously calculated data. 



(f) Definition of set point. 
Figure 4. Continued. 



















* 




* First entry specifies order of numerator, 
second entry specifies order of denominator. 


(i) Definition of parameters for gust approximation. 




< 


RETAIN (15) NORMALIZED GENERALIZED FORCES AS OUTPUTS 


> 




Define node selection array 

(4(D) 


■^ADD (15) SENSORS AT 





(j) Definition of system outputs. 
Figure 4. Continued. 



(1) Define mode selection array. 


Use only if the number of controls is not zero 



(m) Define control selection array. 


Use only if the number of gusts is not zero 



(o) Data card. 
Figure 4. Concluded. 











(1) DAST ARW-1 WIND TUNNEL MODEL, MACH=.9 
(6G12.5) 


(2) MODES, CONTROLS, GUSTS, FREQUENCIES 



10 1 

1 8 





(3) 

GENERALIZED 

MASS 






.21000E-01 

0. 

0. 

0. 

0. 

0. 


0. 

0. 

0. 

0. 

0. 

.44360E- 


0. 

0. 

0. 

0. 

0. 

0. 


0. 

0. 

0. 

0. 

.40220E-01 

0. 


0. 

0. 

0. 

0. 

0. 

0. 


0. 

0. 

0. 

.16960E-01 

0. 

0. 


0. 

0. 

0. 

0. 

0. 

0. 


0. 

0. 

.26920E-01 

0. 

0. 

0. 


0. 

0. 

0. 

0. 

0. 

0. 


0. 

.27170E-01 

0. 

0. 

0. 

0. 


0. 

0. 

0. 

0. 

0. 

0. 


.29440E-01 

0. 

0. 

0. 

0. 

0. 


0. 

0. 

0. 

0. 

0. 

.64510E- 


0. 

0. 

0. 

0. 

0. 

0. 


0. 

0. 

0. 

0. 

.43160E-01 

0. 


0. 

0. 

0. 

0. 

0. 

0. 


0. 

0. 

0. 

.31410E-01 



(A) 

GENERALIZED 

DAMPING 






0. 

0. 

0. 

0. 

0. 

0. 


0. 

0. 

0. 

0. 



(5) 

GENERALIZED 

FREQUENCY 






5.2330 

19.130 

20.910 

25.770 

46.110 

61.230 


79.680 

86.030 

98.090 

118.20 



(6) 

REDUCED FREQUENCIES 






0. 

. 10000 

.30000 

.50000 

.70000 

.90000 


1.3000 

1.8000 





(7) 

GENERALIZED 

FORCE FOR THE, 1 FREQUENCY, REAL 

PART 



25.240 

39.583 

.18119E-01 

-4.9160 

13.585 

3.0988 


4.1912 

-.16004E-01 

-2.9069 

6.4767 

.41256 

-261.26 


-383.03 

-.19168 

31.395 

-64.890 

-31.222 

-72.605 


.20863 

25.085 

-59.666 

-2.7330 

-.20968 

-.32400 


-.15800E-03 

.29615E-01 

-.76497E-01 

-.30163E-01 

-.48386E-01 

.15504EH 


.21652E-01 

-.53256E-01 

-.18149E-02 

-112.69 

-218.92 

-.97377E-I 


24.582 

-109.38 

-26.026 

-3.9417 

.52573E-01 

22.905 


-44.271 

-.20162 

-90.554 

-349.37 

-.14447 

41.350 


-332.36 

-47.701 

-10.520 

.43263E-01 

84.945 

-42.929 


2.1818 

4.8231 

77.844 

.28151E-01 

22.286 

208.21 


-53.838 

-157.73 

.23838 

9.5637 

78.125 

4.1296 


-40.955 

-104.09 

-.61997E-01 

28.705 

88.825 

-58.793 


Figure 5. Example of discrete-frequency data base. Parenthetical numbers (1) to (8) refer to section 6.1 of the 
text. 
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-231.19 

.42353 

-31.645 

33.970 

8.7387 

.67100E' 

.12163 

.89145E-04 

-.24503E-01 

-.22023 

.23544E-01 

.34157 

-.71539E-03 

.19546 

.83258E-02 

-.16810E 

-01 -15.612 

-18.523 

-.14360E-01 

-24.816 

-.37068 

174.85 

131.42 

-.28194E 

-344.77 

-129.31 

7.0539 

9.3481 

21.380 

.86733E 

-8.0774 

49.002 

-9.1440 

7.0522 

.14844E-01 

30.566 

-135.28 

5.3168 

186.44 

-107.84 

-.77173E-01 

-23.090 

-45.771 

-90.593 

-80.140 

.16757 

-100.71 

1.9255 

62.662 

2136.4 

4326.9 

1.9373 

-321.53 

2905.7 

-243.87 

-667.34 

-.17575 

663.29 

647.10 

28.739 

(8) GENERALIZED 

FORCE FOR THE, 1 FREQUENCY, IMAGINARY PART 
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0. 

0. 
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Figure 5. Concluded. 


(1) DAST ARW-1 WIND TUNNEL MODEL, MACH=.9 

(2) NUMBER OF MODES, CONTROLS AND GUST IN THE MODEL 


5 

1 

1 





(6G12.5) 







(3) MODE 

1 

.02366 

1 






.9371 


23.44 

67.66 

16.03 

.4722 

35.23 

33.03 


-22.60 

-.4374 

-4.286 

-14.62 

-11.67 

.4216 


7.995 

-14.86 

-5.188 

-.05001 

4.431 

.5359 


-1.664 





MODE 

1 

.2135 

-53.96 

2 

-233.1 

-106.2 

43.75 

-81.11 

-193.2 

277.8 


-63.34 

5.117 

61.05 

107.7 

-50.92 

-14.13 


69.20 

164.3 

-58.83 

-7.393 

-32.96 

-12.50 

MODE 

1 

.1353 

3 

-7.457 





-15.22 


-92.02 

-77.42 

20.30 

-29.42 

-145.3 

69.22 


-9.049 

3.537 

42.34 

152.2 

4.483 

-15.10 


-48.06 

92.07 

-6.066 

-3.176 

-45.76 

-53.61 

MODE 

1 

.1281 

4 

-3.735 





-10.81 


-52.49 

-75.14 

44.23 

-43.03 

-246.0 

107.2 


1.714 

3.414 

75.90 

231.0 

-84.57 

-45.12 


-208.2 

314.1 

-86.71 

-4.133 

-98.91 

-100.5 

MODE 

1 

.0964 

5 

21.23 





.9128 


-25.74 

-3.885 

-19.63 

6.972 

22.70 

-1.655 


-31.81 

1.707 

13.91 

-130.8 

39.08 

19.39 


164.1 

-99.10 

-7.791 

-7.853 

5.663 

264.2 

CONTROL 

1 

.3336 

1 

-5.296 




1 

61.79 


158.7 

-7.124 

9.744 

-36.74 

-179.4 

-30.73 


15.45 

-7.183 

-29.25 

-21.63 

7.229 I 

-17.91 


-30.03 

13.81 

-19.33 

-31.18 

-89.12 j 

-12.65 


-7.938 





(4) GUST 

1 






0 

1 






2136.4 


4326.9 

-321.53 

2905.7 

-243.87 


(5) GENERALIZED 

MASS 





.21000E-01 

0. 

0. 

0. 

0. 

0. ; 

.44360E-01 

0. 

0. 

0. 

0. 

0. 

.16960E-01 

0. 

0. 

0. 

0. 

0. 

.26920E-01 

.27170E-01 

0. 

0. 

0. 

0. 

0. 

(6) GENERALIZED 

DAMPING 





0. 


0. 

0. 

0. 

0. 


(7) GENERALIZED 

FREQUENCY 





5.2330 


19.130 

25.770 

46.110 

61.230 

79.680 


Figure 6. Example of continuous-frequency data bcise. Parenthetical numbers (1) to (7) refer to section 6.2 of 
the text. 
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Figure 7. Concluded 
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Figure 8. UNAERO structure and hierarchy. 
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